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Abstract: After the recent discovery of a Higgs- like boson particle at the CERN LHC- 
collider, it becomes more necessary than ever to prepare ourselves for identifying its stan- 
dard or non-standard nature. The fundamental parameter Ar, relating the values of the 
gauge boson masses [Jl%, Mz] and the Fermi constant [Gp], is the traditional observable 
encoding high precision information of the Electroweak physics at the quantum level. In 
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this work we present a complete quantitative study of Ar in the framework of the general 
(unconstrained) Two-Higgs-Doublet Model (2HDM). First of all we report on a systematic 
analysis of Ar at the full one loop level in the general 2HDM, which to our knowledge 
was missing in the literature. Thereby we extract a theoretical prediction for the mass 
of the W-boson in this model [M^], taking Mz, a and Gp as experimental inputs. We 
find typical corrections leading to mass shifts <5 M$ ~ 20 - 40 MeV, which help to improve 
the agreement with the experimentally measured value M^ p in a degree no less significant 
than in the MSSM case. In the second part of our study we extend our calculation be- 
yond the mere one-loop order. We devise an effective Lagrangian approach that captures 
the dominant higher order quantum effects on 5p (viz. that part of Ar describing the 
breaking of the approximate SU(2) custodial symmetry) in the limit of large Higgs boson 
self- interactions. This limit constitutes a telltale property of the general 2HDM which is 
unmatched by e.g. the MSSM. Our conclusion is that the Electroweak precision program 
to be conducted at the LHC, and maybe at a future linear collider, can nicely complement 
the direct searches. Should these distinctive loop effects be eventually found they would 
signal a smoking gun hinting at non-standard Higgs physics. 
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1. Introduction 



After the recent announcement of the discovery of a Higgs-like boson candidate at the 
CERN LHC-collider [1-4], we might be closer than ever at unveiling the ultimate archi- 
tecture of the Electroweak Symmetry Breaking (EWSB) mechanism. Since the idea of 
spontaneous symmetry breaking (SSB) was incorporated in the structure of the Standard 
Model (SM), it remained a most pressing and unsettled conundrum in our understanding 
of High Energy Physics. The ideas pioneered by Higgs [5,6], Englert and Brout [7], and 
Guralnik and Kibble [8] crystallized into the present-day paradigm, in which we assume 
the existence of one fundamental scalar SUl(2) doublet whose particle remnant after SSB 
of the gauge symmetry corresponds to the physical Higgs boson of the SM. The quest for 
an experimental confirmation of this picture has hitherto ranked very high in the wish- 
ing list of the experimental efforts first conducted at LEP, and later on at the Tevatron 
and the LHC. The tantalizing Higgs-event candidates reported by the ATLAS and CMS 
collaborations during the last year have apparently been confirmed by the recent analy- 
sis of the new data [1-4], which have led them to conclude (at a ~ 5a confidence level) 
that a new particle has been discovered carrying most of the ingredients to be expected 
from a Higgs-like boson with a mass close to 125 — 126 GeV. The high evidence has been 
gathered by both the ATLAS and CMS collaborations from the detected excesses in the 
diphoton (77) and weak gauge boson (Z°Z°* and WW*) decay modes [1-4]. If confirmed 
by subsequent searches and final identification, this could be the greatest achievement of 
Particle Physics since the discovery of the top quark 17 years ago at Fermilab. And in 
that case it would represent the most impressive and significant success of Particle Physics 
ever, since it would constitute the confirmation of the physical reality of the SSB, i.e. of 
the most subtle and far reaching quantum field theoretical (QFT) structure of the SM. 
Not surprisingly it would at the same time raise many other problems, even outside the 
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strict domain of Particle Physics, such as for example in cosmology through the certified 
existence of the huge electroweak vacuum energy. However, we understand that obtaining 
a consistent overall picture of our world may still take quite some more time. 

In the meanwhile, even after that phenomenal discovery, it is difficult to ascertain the 
nature of the found Higgs-like particle, whether it is the SM Higgs boson or a member of 
an extended Higgs sector whose remaining constituents are yet to be found, or something 
else. There is no reason a priori for scalar particles not to appear within multiple fami- 
lies, as fermions or gauge bosons do. Both theoretically and phenomenologically there are 
many motivations for Higgs physics beyond the SM. For example, the Vacuum Expecta- 
tion Value (VEV) of the Higgs field, v ~ 246 GeV, is well known to be unstable under 
radiative corrections as ultraviolet (UV) contributions to this VEV must be fine-tuned 
to insure its stability at low energies. This naturalness puzzle has constituted a primary 
driving force for the study of extensions of the Standard Model (SM), and of the Higgs 
sector in particular. As a simple, yet very attractive example of such an extension, we 
find the Two-Higgs-Doublet Model (2HDM) [9,10]. Here, by introducing a second SU L {2) 
scalar doublet <3?2 we meet a compelling phenomenological profile [9]- [13]. On a theoretical 
basis, there is indeed strong support for a multi-Higgs doublet structure. To begin with, a 
generic 2HDM furnishes a suitable low-energy description of multifarious EWSB models, 
such as the Minimal Supersymmetric Standard Model (MSSM) [14], Composite Higgs [15] 
and Little Higgs models [16]. Despite its simplicity, it may accommodate a plethora of 
new mechanisms either for spontaneous or explicit CP-violation as well as a rich vacuum 
structure [17]. This sets the ground for novel scenarios in which to address a wide variety 
of unsolved riddles in very different areas, from neutrino mass generation [18] to Elec- 
troweak Baryogenesis [19]. In particular, we cannot exclude more exotic possibilities in 
the framework of Grand Unified Theories, which could also adapt to the properties of the 
purportedly found Higgs boson [20]. 

The simplicity and nevertheless very rich potentiality of the 2HDM qualifies as an 
excellent starting point for a broad class of model-building studies. For instance, if one 
posits additional discrete symmetries in the scalar potential, the generic 2HDM structure 
becomes further restricted and may lead to the so-called Inert Doublet Models (IDM) [21], 
in which at least one of the scalar degrees of freedom becomes a stable WIMP, and thus 
provides a natural candidate for Dark Matter. These models have been portrayed at 
length in the literature, both from the perspective of collider observables (cf. e.g. [22]) 
and astrophysics [23]. A 2HDM structure can also realize a Higgs portal scenario [24]. 
Here, one entertains the possibility of an additional Higgs field in a hidden sector, with no 
couplings to the SM particles - with the exception of the standard Higgs boson. The Higgs 
self- interactions thus constitute the only link (the "portal") between the visible and the 
hidden domains, with a dramatic impact on the expected Higgs boson widths - and so on 
their foreseeable collider signatures [25]. Further extensions of the minimal 2HDM include 
vectorphobic [26] and scale-invariant formulations [27]; or combinations with additional 
gauge bosons [28], fermionic generations [29] or heavy neutrinos [30], among others. 

Numerous studies have scrutinized the prospects for pinning down evidences of 2HDM 
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physics at the LHC, see. e.g. Refs. [31]. In this vein, the potential Higgs-like candidates 
identified by ATLAS and CMS have already endorsed respective analyses on the corre- 
sponding implications for the model [32]. This task is certainly not easy, as long as the 
studies are restricted to direct collider observables. But the situation can improve when 
counting also on the information conveyed by EW precision analyses. In this respect it 
is fair to say that the LHC will soon be able to broaden its capabilities from direct dis- 
covery to precision physics measurements. This is of foremost importance since it may 
provide virtual access to possible new degrees of freedom coupled to the SM, and hence to 
new - or modified - interaction patterns which can induce departures of the EW precision 
observables from the pure SM expectations 1 . 

A most important example of precision observable sensitive to virtual effects from new 
physics is the mass of the charged EW gauge boson [M w ± ]. The current uncertainty in 
its theoretical determination within the SM is estimated to be AM^J ~ 4MeV [38]. On 
the other hand the present world-average of the experimental measurements renders the 
value M^ p = 80.385 ± 0.015 GeV [39], thus carrying an accuracy at a remarkable level 
< 0.02 %, i.e. better than 2 parts in ten thousand. But not less remarkable is the fact that 
the experimental error still gives room for non-negligible non-SM contributions. Indeed, for 
a SM Higgs boson mass of Mr = 125 GeV (cf. the forthcoming discussions in Section |||), 
the corresponding prediction for the mass of the W-boson renders M^ M = 80.363 GeV. 
Despite the current discrepancy | — M^ p | ~ 20 MeV lies to within one-sigma level of 
the experimental measurement, it is as big as five times the estimated theoretical error in 
the SM. Since the latter falls within the reach of the most accurate planned measurements 
of the W-mass in the future there is little doubt that a deviation from the SM of that sort 
should eventually be accessible to observation. This is strongly hinted by the expectation 
that with the upcoming LHC data the uncertainty on the W-boson mass can be narrowed 
down to AM^ P ~ 10 MeV [40], while a high-luminosity linear collider running in a low- 
energy mode at the W + W~ threshold should be able to reduce it even further, namely at 
the A M^ p ~ 6 MeV level or less. The profound impact that a high precision measurement 
of Mw as well as its correlation to the Z-mass, Mz, could have both as a precision test of 
the SM and as a probe of new physics should not be underestimated. This is not new of 
course, what is new is the fact that we are now much closer than ever to exploit this feature 
at the LHC. Let us recall that the M\y — Mz correlation is usually parameterized in terms 
of the quantity Ar, defined as Gf/V% = (g 2 /8 M^)(l + Ar) where Gf is Fermi's constant 



and g is the weak SU(2) gauge coupling- see Sec. 2.3 for a more detailed definition in 
the on-shell scheme. It suffices to say here that g 2 in this scheme is to be replaced by 
e 2 /s 2 N = 47tq/ (l — Mjy/M|), where a is the e.m. fine structure constant. 

The history of Ar and its companion parameter 5p (i.e. that part of Ar which 
parametrizes the breaking of custodial symmetry) is extensive and already quite old 
[33-37]; it has at present more than thirty years. Let us recall that Ar was first computed 
in the SM context in 1980 by A. Sirlin and coworkers [41,42], whereas the 5p parameter 
was defined earlier by Veltmann and collaborators [43-46]. Since then the calculations of 



1 Cf. for instance [33-37] for a comprehensive exposition of the method. 
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these parameters in the SM became improved over the last three decades by important 
QCD and electroweak higher order effects, hence establishing a powerful relation which 
allows to perform accurate predictions of the W-boson mass in high-precision tests of the 
standard model. Not only so, the calculations were soon extended to physics beyond the 
SM, mainly from Supersymmetry (SUSY). In this regard, an exhaustive coverage of Ar 
and bp within the MSSM is available in the literature. The first preliminary calculations 
(including only the so-called oblique contributions) within the MSSM were presented some 
twenty years ago in [47] , and then shortly afterwards at full one-loop level in [48] and [49] 
- see also [50], [51] 2 . Subsequently a lot of refinements which include higher order effects 
have been performed up to the present days, see the comprehensive studies [54], including 
dedicated work e.g. on two-loop effects [55] or flavor violation [56]. 

In spite of the generous literature on the Ar parameter in various contexts of physics 
beyond the SM - mainly focused on the MSSM, as we have just seen - it is a bit sur- 
prising the scarce attention that has been paid to this topic from the viewpoint of the 
general 2HDM, except for some works presented long ago [57-61]. Filling such gap in 
the literature is our main task herewith. To that aim we provide a fully-fledged updated 
analysis of the 2HDM contributions to the parameter Ar at one loop, and use it to ob- 
tain a theoretical one-loop determination of il% in this model, taking Gf, ol and Mz 
as experimental inputs. In addition to that, we extend the pure one- loop analysis in an 
attempt to estimate the maximum impact of some genuine 2HDM higher order corrections, 
namely in the limit where the Higgs boson self-interactions become very large (bordering 
the perturbative unitarity bounds). It is known that these scenarios can strongly modify 
the Higgs/gauge boson couplings, owing to the enhanced quantum effects driven by the 
Higgs self- inter actions. Interestingly, multi-Higgs doublet structures with strongly coupled 
Higgs particles are well motivated from the theory side, in particular in view of strongly- 
interacting realizations of the EWSB. Our main focus in this second part of the paper will 
be to determine whether such augmented Higgs self-couplings may be able to stamp any 
sensible fingerprint on the EW precision observables under scope. 

The paper is organized as follows. In Section ^ we succinctly review the 2HDM setup 
and relevant constraints; we also consider a preview of the EW precision observables to be 
examined in detail thereafter, and set our definitions and notation. Section || is devoted to 
present the results of our detailed analysis of the predictions for Ar and M\y at full one 
loop level in the 2HDM. In Section ^| we explore the 2HDM region beyond one- loop level 
by introducing an effective Lagrangian approach that enables us to cope with the leading 
higher order effects on 5p in the large Higgs self-coupling limit. Conclusions and closing 
remarks are finally delivered in Section ||. Additional analytical details of the calculation 
are quoted in the Appendix. 

2 Let us point out that to the best of our knowledge the oldest full one- loop MSSM calculation of 
the electroweak gauge boson masses existing in the literature was provided quite earlier in references [52] 
and [53]. Although it was presented in a renormalization framework slightly different from the usual one, 
it was later adapted to the the standard on-shell scheme and this resulted in the first full one-loop MSSM 
calculation of Ar reported in the literature [48] , followed shortly after by the similar analysis of [49] . 
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2. Theoretical setup 



2.1 The general Two-Higgs-Doublet-Model in a nutshell 

The 2HDM [9] canonically extends the SM Higgs sector with a second SUl(2) doublet of 
weak hypercharge Y = +1, so that it contains 4 complex scalar fields. The most general 
form of a gauge invariant, renormalizable, CP-conserving potential that one may construct 
out of two doublets $j (i = 1, 2) can be casted as follows 



where the self-couplings Aj may be rewritten in terms of the masses of the physical Higgs 
particles (M h o, M H o, M A o, M H ±); tan/3 = i^/ui (the ratio of the two VEV's (0°) giving 
masses to the up- and down- like quarks); the mixing angle a between the two CP-even 
states; and, last but not least, the self-coupling A5, which cannot be absorbed in the 
previous quantities. Therefore we end up with a 7-parameter set: (M h o, M H o, M A o, M H ±, 
sina, tan/3, A5). An additional discrete symmetry <l>j — > (— l)*<3?j(i = 1,2) - which is 
exact up to soft-breaking terms of dimension 2 - is canonically assumed as a warrant of 
Flavor-Changing Neutral-Current (FCNC) suppression 3 . Alternative constructions with 
no explicit Z 2 symmetry are described e.g. in [62] and references therein. 

As for the Yukawa sector involving the Higgs/quark interactions, the absence of tree- 
level flavor changing neutral currents (FCNC) leads to two main canonical realizations: 1) 
type-I 2HDM, in which just one Higgs doublet couples to all quarks, whereas the other 
doublet does not; 2) type-II 2HDM, where one doublet couples only to down-like quarks 
and the other doublet just to up-like quarks. Other flavor structures are also conceivable 
and have indeed attracted a growing attention in the recent years [63]. That said, we will 
see that the evaluation of Ar is barely influenced by the particular form of the Yukawa 
couplings after we impose the various phenomenological and theoretical restrictions. For 
a more detailed exposition within our notation, including a detailed list of Higgs boson 
couplings to fermions and bosons, and a discussion of the on-shell renormalization of the 
unconstrained 2HDM Higgs sector, see Ref. [64]. 

It is also worth recalling that the Higgs sector of the MSSM corresponds to a partic- 
ular (constrained) realization of the general two-Higgs doublet structure. The underlying 
SUSY invariance restricts the form of the potential (|2.1| ) in a way that has far-reaching phe- 
nomenological implications. While in the MSSM the Higgs self- interactions, and thereby 
also the Higgs mass spectrum, are dictated by the EW gauge couplings, in the 2HDM we 
no longer have such a dynamical restriction a priori. This implies that the triple, as well as 
the quartic, 2HDM Higgs self-interactions are fundamentally unconstrained - the influence 
of the latter being comparatively milder. In practice, they can be boosted as much as per- 

3 Let us note in passing that such a symmetry is automatically preserved in the MSSM. 
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mitted by perturbative unitarity and vacuum stability, giving rise to trademark signatures 
which are completely foreign to the MSSM 4 . 

A rather extensive set of bounds restricts the regions of the 2HDM parameter space 
with potential significance to Phenomenology Dedicated accounts on these topics can be 
found e.g. in Refs. [62,69-71], and in particular also in Section II of Ref. [64], which very 
closely follows the notation and conventions employed herewith. It goes without saying 
that the recent identification of a ~ 125 GeV SM-like Higgs boson candidate by ATLAS 
and CMS raises a number of very significant implications, which we must take into account 
for realistic studies. A comprehensive updated analysis of the 2HDM parameter space 
constraints in the light of these novel results is not yet available. Nevertheless, we can rely 
on model-independent approaches [72] that spell out the general conditions to be satisfied 
by phenomenologically viable extensions of the standard Higgs sector. 

2.2 Phenomenological restrictions 

The basic restrictions to be imposed on the 2HDM parameter space can be outlined as 
follows: i) first, we need to account for an upper limit on the loop-induced breaking of the 
(approximate) SU(2) custodial symmetry, which one can precisely trade through Sp, that 
is, one of the EW precision quantities at stake in this paper. Its numerical value should lie 
below I^hdmI < 0(1(T 3 ) [39], if we allow up to 3a deviations from the current best-fit 
value. In practice, this translates into restrictions on the mass splitting among the different 
Higgs bosons. Tight requirements ensue also from the radiative B-meson decay B(b — > 57) 
and the - B° d mixing. Generically, the former process sets a lower bound on the charged 
Higgs mass of M H ± > 300 GeV for tan /3 > 1 [62] in the case of type-II 2HDM, while the 
latter strongly disfavors the tan/3 < 1 regions (for both type-I and type-II) and tends to 
enforce tan /3 to be roughly above 1.5 — 2 if the charged Higgs boson is kept relatively light 
(viz. M H ± ~ 100— 150 GeV, which is allowed for type-I models). Perturbative unitarity, as 
well as vacuum stability, impose as well very severe limitations. These translate into wide 
excluded areas across the tan/3 — A5 plane. Unitarity places an upper limit of |As| ~ O(10) 
for Higgs boson masses of few hundred GeV, and vacuum stability excludes the A5 > 
region up to a very narrow band [64]. In short, in order to satisfy these restrictions we are 
confined to regions where tan/3 ~ 1 and [As) ~ 5 — 10 (A5 < 0). Additionally, any chosen 
Higgs mass spectrum ought to satisfy all the current limits from direct searches at LEP, 
Tevatron and LHC. 

In practice, all these constraints are systematically included in our calculation by 
combining the latest version of the public codes 2HDMCALC-1.1 [73], SuperISO-3.1 [70] 
and HiggsBounds-2.2 [74], hand in hand with several alternative and/or complementary 
in-house routines. Last but not least, we must deal with the phenomenological implications 

4 In contrast, the core of the enhancement capabilities of the MSSM Lagrangian resides in the richer 
pattern of Yukawa-like couplings between the Higgs bosons and the quarks, as well as between quarks, 
squarks and charginos/neutralinos. Their implications for collider and EW precision physics have been 
object of dedicated attention in the past for a plethora of varied processes, see e.g. [65]. For reviews on the 
subject, see e.g. [9,66-68]. 
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of the 5a Higgs boson candidate recently unveiled at the LHC [1—4]. A most natural 
choice when embedding the current experimental picture into a concrete realization of 
the 2HDM is to identify the lightest neutral CP-even Higgs boson [h°] with the ~ 125 
GeV resonance. Under this assumption, the Higgs couplings to the gauge bosons become 
severely constrained, as the current data show no substantial departure with respect to 
the SM-like decay patterns. Consequently, we are left with very tight restrictions on the 
trigonometric factors tan/3 and sin a. It follows that we are essentially restrained to the 
so-called decoupling regime a = (3 — ir/2 (with a < 0) or, equivalently, a = (3 + ir/2 (with 
a > 0) - both cases featuring g^o vv ~ sin 2 (/3 — a) ~ 1. If we instead identify the ~ 125 
GeV resonance with the heavy neutral CP-even Higgs boson [H°], the observed decay 
rates into gauge bosons enforce a ~ (3 since then, again, the corresponding coupling is 
maximized: g%o vv ~ cos 2 (/3 — a) ~ 1. Independently, these conditions also disallow some 
specific choices of sin a. For instance, sin a ~ cannot be realized within a type-II 2HDM. 
In the decoupling limit, this choice would imply (3 ~ — 7r/2, thus rendering an unduly 
enhanced h°bb interaction, ~ | sin a/ cos (3\ » 1, incompatible with the fermionic modes of 
the current Higgs candidate observations. Conversely, for a = (3 the choice a ~ would 
spoil the perturbativity of the Higgs/top Yukawa coupling. By similar arguments, one 
can prove that a ~ 7r/2 is not permitted within type-II realizations of the 2HDM. Notice, 
however, that a = (resp. a = vr/2), namely the fermiophobic limit for H° (resp. h°) is 
still viable within type-I models owing to the different Higgs/quark interaction strengths 
~ sin a/ sin/3 (resp. ~ cos a/ sin/3). 

Further constraints can be imposed over tan/3. For example, we must comply with 
the aforementioned B-physics limits on both the low (tan/3 < 1) and high (tan/3 > 10) 
tan/3-ranges. Also the non-observation of enhanced h — > bb decays rules out the large 
tan/3 regimes of a type-II 2HDM. In addition, recent studies (cf. Ref. [72]) have concluded 
that the fermionic decay signatures of the ~ 125 GeV Higgs-like resonance, even if fully 
compatible with a SM-like Yukawa sector, exhibit a mild statistical tilt towards slightly 
enhanced (resp. suppressed) Higgs/top (resp. Higgs/bottom) couplings. These scenarios 
would be realized for tan (3 < 1. Let us recall, finally, that moderate choices of tan /3 ~ 0(1) 
are particularly appealing and well motivated for the purposes of our analysis, as they 
enable to maximize the triple (3H) self-couplings through a relatively large value of the 
parameter | A5 1 - therefore by resorting only to the intrinsic structures of the 2HDM Higgs 
potential. 

A variety of processes can probe the potentially enhanced 2HDM Higgs boson self- 
interactions, and have indeed been intensively analysed over the past years - mostly in 
the context of linear colliders. Available studies include, on the one hand, the tree-level 
production of triple Higgs-boson final states [75]; the double Higgs-strahlung channels 
hhZ° [76]; and the inclusive Higgs-pair production via gauge-boson fusion [77]. In the 
same vein, also the 77 mode of a linac has been explored, in particular the loop-induced 
production of a single neutral Higgs boson [78] and of a Higgs boson pair [79]. In all the 
above mentioned cases, promising signatures were pinpointed, which could be revealing of 
an unconstrained multi-Higgs doublet pattern. Similar genuine 2HDM effects might also 
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manifest as large radiative corrections to a number of Higgs production channels. One loop 
studies of pairwise Higgs boson final states were first addressed for charged Higgs bosons 
e + e~ H + H~ [80] and later on carried to completion by a full-fledged study of the neutral 
Higgs sector e + e _ — > h°A°, H°A°, including also the more traditional Higgs-strahlung events 
e + e~ — > h Z°,H°Z [64,81]. All these studies reveal the possible existence of: i) sizable 
Higgs boson production rates, typically in the ballpark of O(10 — 100) fb for ^/s = 500 
GeV at a future linear collider; ii) sizable quantum effects, up to 5 ~ ±50%; and iii) a 
very characteristic complementarity of the dominant Higgs production modes at different 
center-of-mass energies - these properties being correlated, once more, to significant Higgs 
boson self-interactions, and so to a non-supersymmetric multi-Higgs doublet structure. 




2.3 Electroweak Precision Quantities from muon decay 

The relation between the EW gauge boson masses (Mw, Mz) in terms of the Fermi constant 
(Gf) and the fine structure constant (a) is an essential tool for testing the quantum effects 
within the SM as well as to place bounds on its manifold conceivable extensions. Such 
relation can be derived in terms of the muon lifetime r^, whose decay rate is precisely 
defined by the Fermi constant, Gf, via the expression [33-36] 



1 + A Q ed), (2.2) 



where F(x) = 1 — 8x — 12x 2 lnx + 8x 3 — x 4 . Following the standard conventions in the 
literature, the above defining equation for Gf includes the expression Aqed, i-e., the finite 
QED contribution obtained within the Fermi Model, which is known up to two-loop order. 
Calculating the muon lifetime within the SM at the quantum level and comparing with 
Q2.2| ) yields the relation: 

in which 

Ar ^^wW +Ar [ve rt ,box]_ {2A) 

These expressions define the quantity Ar in a precise way. Let us notice that £w(^ 2 ) is the 
on-shell renormalized self-energy of the W-boson; it accounts for the universal ("oblique") 
part of the electroweak radiative corrections to the muon decay. The non-universal (i.e. 
process-dependent) corrections - which stem from the vertex and box contributions to the 
muon decay - are encoded in the subleading term Ar' vert,bo>! L The explicit expression for 
Ar consists of a combination of loop diagrams and counter terms. For the renormalization 
details, see e.g. [33-36]. Here we will only remind the reader of some basic facts which 
can be helpful to contextualize our 2HDM computation. To start with, let us write down 
the explicit structure of Ar after renormalization: 
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(2.5) 



st^ \ M 2 Ml Ml % Ml 



We stress that this expression is finite because its original definition fl2.4|) depends only 
on the on-shell renormalized self-energy of the W-boson and the remainder Ar' vert ' box ' - 
which is also finite by virtue of the Ward-Takahashi identities of the EW theory. In the 
previous expression II 7 (A; 2 ) = d'E^ f (k 2 )/dk 2 is the photon vacuum polarization, and we 
are using the notation si = 1 — M^/M|, and cl = M(y/M|. Furthermore, since the 
renormalization is in the on-shell scheme the weak gauge boson mass counter terms read: 
5M 2 = 3?e£y(M 2 ), with V = W ± ,Z°. By £y(g 2 ) = T^{q 2 ) we will hereafter denote the 
(transverse parts of the) unrenormalized gauge boson self-energies, which we conventionally 
extract from the corresponding vacuum polarization tensor: 

IW) = gT ^W) + <t <f A (Q 2 )- (2-6) 

After introducing the renormalized photon vacuum polarization II 7 (A; 2 ) = II 7 (A; 2 ) — II 7 (0) 
it is convenient to consider its fermionic part at k 2 = M§, i.e. II 7 rm (M|). As it is well- 
known, for the light fermions (quarks and leptons, excluding the top quark) the quantity 

Aa = -Ken^ rm (M|) (2.7) 

is independent of the EW part of the SM and goes into a finite renormalization of the 
QED fine structure constant: a —> a(l + Aa). This can be resummed according to the 
renormalization group to provide the value of a at the scale of the Z-mass: a(M|) = 
a/(l — Aa) ~ 1/128 which is ~ 6% larger than its low energy (Thomson limit) value 
a ~ 1/137. Such resummation takes into account all the leading logarithms of the type 
a n ln n (Mz/mi) which enter the renormalization of a from the leptonic sector I = e, /i,r. 
The light quark contribution, instead, is computed more precisely via a dispersion relation 
from the experimental hadronic data collected in low energy e + e~ scattering. Finally, the 



top quark gives a negligible (decoupling-like, ~ M^/m 2 ) contribution to (2/7). 

Let us now emphasize one more type of effect which will have special relevance for our 
analysis of the 2HDM contributions to Ar. We are referring to the so-called bp parameter 
[43-46]. By inspecting fl2.5| ), such contribution comes from the term 

dp (2.8) 



5M£_6M1 s z (Q) s w (Q) 
M| Ml M 2 Ml 



The 5p effect is finite for the matter fermions of the SM (for each doublet separately) 
and provides a very important (non-decoupling) contribution when large mass splitting are 
present within a given fermion family. Its main source comes of course from the top quark, 
or to be more precise from the top quark and bottom quark doublet. Even though the 
bottom quark contributes a negligible finite correction, it is essential to make the quark 
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doublet contribution to Sp perfectly finite on its own. The corresponding impact on Ar 
is not just — (c^/s^y) 5p l as it actually contains additional terms (encapsulated in the the 
so-called "remainder" Ar rem , cf. Eq. ( 2.11J) and the discussion further down) which are 
numerically significant: 

Ar top = 



V2G F M^ 



16vr 2 



rrif 4 2 c w 7 

3 4^ +2 l4;"3j ln ^; + 3 lnCw + 4;-9 




(2.9) 



For a more physical interpretation, let us remind the reader that 5p stands for the possible 
deviations of the value of the Fermi constant in neutral current processes (Gp C ) - typi- 
cally induced by neutrino interactions - from the corresponding Fermi constant in charged 
processes (Gp) such as muon decay. Both in the SM and in the 2HDM we have 

'=*5r = iSr 1+ *- <2 ' 10) 

In the absence of weak hypercharge interaction g' — > (Mz — > M\v, S\v ~~ ^ 0) we WOUi d 
have (5/9 = and the Fermi constants in both kind of processes would have equal strength. 
In the SM, as well as in the 2HDM, the only sources of Sp come from quantum effects. 
These deviations are bound to satisfy \Sp\ < 0(1O~ 3 ) [39]. But this bound still gives a 
substantial margin for physics beyond the SM, as we shall see. 

After clarifying the physical significance of the main terms in ( p.5| ) we may now write 
down the general structure of Ar in the following traditional form [33-36]: 

c 2 

Ar = Aa - 5p + Ar rom = Aa + Ar^ + Ar rem , (2.11) 

where the leading contributions Aa and 5p have been defined above, and we have intro- 
duced Ar^ p ] = — (cyy/syf)5p. The so-called "remainder" piece Ar rem collects the remaining 
effects, which in the SM entail subleading (which should not be taken as synonymous of 
negligible) contributions . For example, while Aa — 0.06, we have Ar rem ~ 0.01. For 
comparison, the top quark gives a contribution to ( pUD of around Ar^]( top ) ~ -0.03, 
or —0.04 if using the more accurate expression fl2.9|) . Although Ar rem is smaller than the 
leading terms, all these contributions are in fact quite significant and produce an impor- 
tant numerical shift of the W-mass, lowering its zeroth order value (M^p e ~ 80.9 GeV) by 
about 0.8%, i.e. more than 0.6 GeV. In particular, Ar rem renders a non-negligible contri- 
bution of ~ 160 MeV to that total (see below). The additional effects that might come 
from physics beyond the SM are generally much smaller but the current high precision 
electroweak physics has improved sufficiently so as to make possible to detect shifts of Ar 
at the level of less than one per mil. At the end of the next subsection we shall further 
illustrate this point. 

At present, the calculation of Ar in the SM is complete up to two loops [55,83-86] and 
includes also the leading three and four-loop pieces [87]. As we have mentioned, the bulk 



5 Beyond the one- loop order, resummations of the leading one-loop contributions have been derived, see 
e.g. [82], and further contrasted to exact higher-order calculations [55,83]. 
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of the contributions to Ar is linked to the renormalization of the fine structure constant. 
The next-to-leading source of contributions comes from Sp, which in the SM is finite and 
dominated by the aforementioned O(mf) terms from the top-quark loops. In the SM, 
however, the Higgs contribution to 5p is neither gauge invariant nor UV-finite, only the 
sum with the remaining bosonic part is finite and gauge independent. Numerically it is 
not very relevant as compared to the fermionic contribution to Ar and it increases only 
logarithmically with the Higgs boson mass. This feature reflects the so-called screening 
behavior of the SM Higgs boson [44], a property which does not generally hold in extended 
Higgs sectors, such as e.g. in the unconstrained 2HDM. As this issue is important for 
our considerations, let us quote explicitly the leading effect of the SM Higgs boson in the 
limit of large Mjj- The contribution being not finite, it also depends on the choice of 
the dimensional regularization scale p used in the 't Hooft-Veltman procedure (see the 
Appendix for more details). The splitting of terms in the full bosonic contribution is 
therefore somewhat arbitrary, but it is natural to set that scale at the EW value p = Mw- 
Finally, omitting the UV-divergent pieces that cancel with the remaining bosonic terms 
one arrives at 



SaP ~ =— ^ln— } + ... (2.12 



and 

H V2G F M^ ii ( m 2 h 5 



U>7T- 3 [ M^r 6 



Ar*~ — , - ^ ln^f -- n +••• (2-13) 



As we can see the dominant term —(c^/s^)5p H is corrected by non-negligible additional 
finite parts from Ar rem having the same structure. The neat contribution of Eq. ( 2,13| ) is 



subsumed again in Ar rem . Let us indeed note that numerically Ar H is rather irrelevant 
as compared to, say, the top quark contribution and the overall Ar value (~ 0.04) within 
the SM; we find Ar H = 0(1O~ 3 ) for M H = 200 GeV and 0(1O~ 4 ) for M H = 125 GeV. 
In spite of this meager yield, we wish to stress that the particular 5p H piece in Eq. ( 2.12; ) 



is formally very important because it measures the departure from custodial symmetry, 
which is that global SU(2) symmetry of the Higgs SM Lagrangian which is only broken 
by the weak hypercharge interaction g' = gsw/cw- The "custodial symmetry limit" thus 
corresponds to g' — > (Mz — > M\v, s 2 ^ ~~ ^ 0)- I n this limit the gauge bosons and 
Z form a degenerate triplet of a global SU(2)r symmetry (custodial symmetry [46]). We 
may indeed confirm from the above expressions that the 5p H effect takes on the form 
5p H ~ (— 3(/ 2 /167r 2 ) In Mjj/Myy and hence it vanishes in the custodial symmetry limit 
g' 2 = g 2 s 2 N /c 2 N — > 0, as expected. 

A natural question to ask is if there are significant (non-screening) custodial-breaking 
effects from Higgs physics beyond the SM, i.e. effects not just growing logarithmically with 
the Higgs boson masses but as powers of the masses themselves. Let us briefly mention the 
Higgs sector of the MSSM. Although it features a type II 2HDM, it is of a constrained nature 
owing to the underlying SUSY [9,68]. As a result there are no conspicuous non-screening 
effects and in this sense the situation does not differ significantly from the SM. There 
are notwithstanding alternative contributions to 5p in the MSSM which come from the 
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stop/sbottom-mediated quantum effects and are potentially very relevant [11,12,47,48,50]. 
Finally, in the general (unconstrained) 2HDM case the situation changes dramatically. The 
corresponding bp contributions have been studied at one-loop [9,57-61,88] and have been 
employed in studies of combined 2HDM parameter space constraints [71]. However, to 
the best of our knowledge no systematic study of the full one-loop contributions to Ar 
from the 2HDM is available in the literature. We believe that in the light of the present 
experimental situation at the LHC it is highly desirable to cover this gap. 

2.4 Ar, bp and Mw in the 2HDM. Preliminary considerations 

It is important to emphasize that Ar is sensitive (through quantum effects) to all the SM 
parameters (couplings and masses) as well as to the full list of parameters involved in any 
possible extension of the SM. In the relevant 2HDM case under consideration, 



Ar = Ar(e, M\y, M%, uif] Mi, sin a, tan f3, A5) 



(2.14) 



where m f and Mj = M h o , M H o , M40 , M H ± are respectively the masses of the fermions and 
of the Higgs bosons. Finally, A5 is the parameter of the Higgs potential Q2.1| ) which cannot 
be related to a physical mass of the 2HDM. That parameter enters Ar beyond one loop 
only, but we will see that it can still furnish significant quantum effects. In general these 
effects will enter Eq. ( 2.11| ) through bp and Ar rem - Let us clarify that the contributions 
to Ar from the charged Higgs bosons come only from bp and Ar rem in Eq. (2.11), as the 
charged Higgs boson loop effects on Aa leave no finite yield to the photon self-energy nor 
to the photon-Z mixing, as expected from gauge invariance. 

In this work we revisit the calculation of bp in the general 2HDM and shall consider 
also the full one-loop calculation of Ar, i.e. the complete quantity (|2.5|) or (2.11). As 



advertised, we will subsequently include higher order effects on bp related with the self- 
couplings of the Higgs bosons, mainly governed by the A5 parameter. These effects are 
non-existent in the MSSM owing to the pure gauge nature of the self-couplings in the 
tree-level MSSM Higgs potential [9]. 



The dominant part of the 2HDM corrections comes from the one-loop contributions to 
bp mediated by the various 2HDM Higgs bosons [57-61]. As in the case of the SM we 
perform the calculation in the Feynman gauge (cf. Fig. |lj) and we call the result bp2HDM- 
It amounts to a compact and finite expression whose full form is provided in the Appendix, 
see Eq. fl6,6|). Here we single out only the portion of 5/J2HDM that depends on the Higgs 
boson mass splittings, which may be called <5p2HDM- It can be casted as follows: 



$P*2 
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Figure 1: Set of Feynman diagrams describing the pure 2HDM one- loop contributions to the 
and Z° self-energies in the 't Hooft-Feynman gauge. All Feynman diagrams in this paper we have 
generated by means of Feyn Arts [89] . 



From this expression, and by comparison with Eq. ( 2.121 ), it is clear that in the uncon- 



strained 2HDM the contributions to bp do not follow the screening theorem of the SM 
Higgs boson. Indeed, from the quadratic dependence on the various Higgs boson masses 
outside the logarithms and the fact that arbitrary mass splittings between these Higgs 
bosons are possible, one can expect that they could easily overshoot the limits on bp 
within the 2HDM, namely the 3er bound \bp\ < O(10~ 3 ) obtained from the EW precision 
fits [39]. However, we note that if M^o — > Mjj± then <5p2HDM — > 0. Hence if the mass 
splitting between M A o and M H ± is not too large ^p^hdm can ^ e kept un der control. This 
is also true for the full <5/52HDM) if at the same time we are not far away from the decoupling 
limit j3 — a = ir/2 - where the lightest CP-even state h° behaves SM-like (see Appendix 
for more details). 

Let us now provide some discussion on the general strategy we will follow to estimate 
Ar in the unconstrained 2HDM at one- loop and beyond. More detailed considerations 
will be made in sections 3 and 4, together with the full quantitative analysis of course. 
Let us assume that Ar in the defining equation ( |2.3j ) collects the expanded form of the 
various contributions up to some order of perturbation theory where they are presumably 
computed. We can split Ar2HDM as follows: 

Ar, mM - Ar [1](RG) + Ar (boson) + Ar (Ykw " QCD) (2 16) 

"'2HDM 2HDM ' 2HDM ' 2HDM ■ V^- iL V 

This structure assumes that all of the leading (and some next-to-leading) effects on Ar 
are involved within the 2HDM. There are various subclasses of contributions in it that 
will deserve some particular comment. At the moment we note that the first term on the 
r.h.s. of ( |2.16| ) is the improved one-loop correction, namely the total one- loop correction 
plus the renormalization group (RG) resummed effects of Aa, which go into the important 
renormalization of the e.m. fine structure constant; the second term represents the complete 
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Figure 2: Small representative sample of Feynman diagrams describing potentially significant 2- 
loop contributions to Ar within the 2HDM in the 't Hooft-Fcynman gauge. The complete list of 
2-loop diagrams of this class is much more extensive, and the full list of 2-loop graphs of all classes 
within the 2HDM (involving bosons and/or fermions) is huge. However, even with the displayed 
sample we can immediately appreciate the presence of trilinear Higgs boson couplings and their 
foreseeable significance. In the limit of large A5, all these diagrams are of C(a ew A|) and constitute 
the (gauge invariant and finite) leading 2-loop contribution to Ar. 



2-loop bosonic part, i.e. the pure EW contribution from gauge bosons, Goldstone bosons 
and Higgs bosons in all possible combinations; finally, the third term stands for the joint 
Yukawa-coupling and QCD effects to 0(a ew a s ), including some resumed contributions. 
Needless to say all these terms are rather complicated but many of them exhibit a structure 
very similar to that of the known SM case [90]. There is, however, a subset of higher order 
diagrams in the 2HDM which provide a qualitatively (and maybe also a quantitatively 
important) new contribution. The latter is buried in the Ar^^^ piece of ( fTp and for 
this reason it proves convenient to further split it as follows: 

A (boson) _ .(boson*) . (A 5 ) ,„ , y x 

2HDM 2HDM ' 2HDM " \*- L <) 

The special term here is the second one on the r.h.s. of the above expression, Ar^^j^, 
whose relevance will soon become apparent. We start by explaining the meaning of the first 
term, Ar^j^j^, which stands in part for the ordinary gauge boson and Golstone boson 
diagrams up to two loops. Among these diagrams we find the 2-loop pieces contributing to 
0(a!g W ), which are not particularly important from the numerical point of view, although 
they are technically difficult to account for. Furthermore, A^hdm*'' involves as well the 
2-loop diagrams containing Higgs bosons only, and also those graphs mixing gauge bosons 
and Higgs bosons; the corresponding amplitudes can be constructed from the one-loop 
diagrams of Fig. [I] after bridging them with another Higgs internal line (propagator). See 
Fig. ^ for a detailed sample of 2-loop diagrams constructed in this way. The graphs in that 
figure are proportional to the product of two trilinear Higgs boson self-couplings, \hhh-, 
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and two ordinary gauge couplings, i.e. the corresponding amplitudes are of 0(a ew Xf lhh ). 
There are other 2-loop graphs constructed from the aforementioned bridging procedure 
that involve three gauge couplings and one trilinear Higgs coupling, and hence are of 
0(g 3 Xhhh)- We can expect that some of these amplitudes, especially those involving two 
trilinear couplings, can be very important from the quantitative point of view because 
they can be enhanced even after preserving all known bounds on perturbative unitarity, 
custodial symmetry and vacuum stability. For a full display of the detailed structure of 
the trilinear Xhhh self-couplings in the general 2HDM, see e.g. Ref. [64] and particularly 
its Table II. As an example we quote the coupling of three light CP-even Higgs bosons: 



A h o h o h o = -— — [M^ (2cos(a + /?) + sin 2a sin(/3 - a)) 

ZMw sin zp s\v 

— cos(a + p J cos (p — a) 7, 



(2.18) 



These structures can be enhanced, in principle, on three accounts: a) assuming large 
tan/3 > 1 or small tan/3 C 1, b) arranging for large mass splittings among the Higgs 
bosons, and c) pushing the value of the characteristic A5 coupling (the free parameter which 
is unrelated to masses) in the general 2HDM Higgs potential ( |2.1j) . Of these possibilities 
in practice only the last one can be used after enforcing all the known theoretical and 



phenomenological bounds described in Sec. 2.2. 

Let us note in passing that the A HHH contributions are not present in the SM at 
two loops (with H standing here for the SM Higgs boson) because the Z gauge boson 
cannot couple to HH in the SM and hence a double HHH vertex at two loops cannot be 
formed by the aforesaid bridging procedure. Only a single vertex HHH can appear at this 
order in the SM, thus giving rise to 0((7 3 Ahhh) effects from the Higgs self-interaction in 
combination with pure gauge boson couplings. But even in this case Ahhh has no special 
enhancing property in the SM beyond the artificial raising of the Higgs boson mass, recall 
that Ahhh ~ gM^/My?. For a relatively light Higgs boson, Mh > 125 GeV the trilinear 
Higgs boson self-coupling in the SM is not very substantial. In particular, for a 125 GeV 
Higgs boson we get 



3eM£ 



Ahhh(M h = 125 GeV) - ., 

z A2\v S 



W 



-188 GeV. (2.19) 

M„=125 GeV 



To certify the meager impact of this numerical result notice that it yields roughly 2% of 
the Lee-Quick-Thacker bound only [91], which as we know represents the upper bound on 
the SM Higgs boson self-interaction ensuing from the perturbative unitarity requirement. 
In view of the fact that the presumed value of Mh is not much higher than Myy, the 
mentioned effects can be finally subsumed into the fairly irrelevant pure bosonic 0(ag W ) 
part in the context of the SM [84,85]. 



Let us now finally address the precise origin of the Ar^pj^ piece in Eq. ( 2.17 ). Such 



contribution appears when we select the A| tag in the subclass of all the 2-loop diagrams 
containing two trilinear Higgs boson couplings. Recall that the A5 term is present in all the 
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Higgs boson self-couplings Xhhh- (At this important point we refer the reader once more 



to Table II of Ref. [64]). In particular, see the last term on the r.h.s. of Eq. ( 2.18j ). It was 
precisely in order to single out such peculiar and potentially relevant 2-loop A| proportional 
effects in the 2HDM computation of Ar that we have introduced the last term of Eq. ( 2.17] ) . 



In our approach we wish to explicitly detach this term from the 2-loop bosonic contribution. 
To do this in practice we split all the trilinear Higgs boson self couplings of the 2HDM as 
follows 

Xhhh = X Hhh + Ch\$ (2.20) 

where X* hhh is the part of the coupling not containing A5, and Ch is a coefficient which 
can be easily identified from the list of trilinear couplings in Table II of [64]. Let us now 
emphasize that the part of the amplitudes which is proportional to X* hhh is also incorporated 
into A^2HDM*^ in Eq. ( 2.17] ) whereas the part which carries the A§ coefficient constitutes 



our precise definition of the piece Ar^DM in ([2.17 )■ In practice, we shall project the 



0(A§) component for all of the Higgs boson self-couplings from the generic structure ( 2.2C| ), 
meaning that, at the end of the day, Ar^ 5 -^ will be built upon O(X^) contributions; these 
are the leading (Higgs-mediated) effects beyond 1-loop in the limit of large Higgs boson 
self-interactions - for more details cf. our discussion in Section ||. 

It is important to note that the formal separation of these two kind of contributions 
is well-defined, and it is worthwhile performing it because the As-part is potentially the 
most significant one from the phenomenologically point of view. We note that even though 
both X* hhh and Ch could also be enhanced at large or small tan (3, we do not contemplate 
this possibility in view of the perturbative unitarity constraints, alongside with the phe- 
nomenological preference for tan f3 ~ 0(1) values. As a result the entire enhancing power 
of the Xhhh couplings resides exclusively in the A5 part, whereas the rest of the 2-loop 
bosonic yield is subsumed in Ar^^*^ , altogether quite irrelevant in practice as in the SM 
case [84,85]. 

The sum of all the 2-loop bosonic diagrams involved in the computation of Ar^^j^ 



in Eq. ( 2.1 7| ) is well defined, namely it is gauge invariant and finite. Most important for 



our purposes is that we have checked explicitly the gauge invariance and finiteness of the 
overall piece Ar^ 5 ^ constructed from the above procedure (see section 4 for more details). 
This feature could be expected because the A5 coupling does not appear at one-loop in this 
calculation and hence that parameter is not needed for the renormalization of the gauge 
boson sector at 2-loops. As a result we can exploit this property and test the impact of the 
2-loop diagrams involving one or two trilinear Higgs boson couplings in the limit of large A5 . 
Of course "large" means, in this context, as large as allowed by the perturbative unitarity 
constraints, which are well-known in the literature [92] - see e.g. [64] for a discussion in our 
framework - and we certainly impose them as an important restriction in our calculation. 
From the above considerations the following strategy is suggested to obtain an esti- 



mate of Ar in the general 2HDM, Eq. (|2.16[) . First, the pure bosonic higher order effects 



independent of A5 are as tiny in the 2HDM as they are in the SM; and, second, in the 
regime tan/3 = 0(1) we should have in the 2HDM at most the same Yukawa couplings 
and QCD loop corrections as in the SM [55,83-85]. We express these two statements in a 
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nutshell as follows: 



A (boson*) _ a (boson) _ » (Ykw-QCD) _ » (Ykw-QCD) ,„ ^ 

2HDM — ^'SM ' 2HDM — ^' SM ' 

^(boson*) 



where the bosonic part At^hdm is of course the same quantity that we defined in ( 2.17 ) 
through the splitting (2.20), and Arg^ son ^ is the known 2-loop bosonic SM result. A 



practical recipe for a reasonable estimate of Ar in the general 2HDM should therefore be 
the following effective quantity: 

Ar^ DM ^ Ar SM + S(Ar^ DM ) + Argf DM . (2.22) 

The first term, ArgM, stands for the full set of SM contributions known to date, hence 
with a structure completely similar to ( 2,16| ) but within the SM: 



ArsM = ArWf G > + Ar^ + ArgT"^ . (2.23) 

It includes different sorts of one loop and higher order effects of various kinds [55,83-85], 
as well as some resummed contributions - the most important one being of course the one 



affecting Aa. The overall SM contribution (2.23) is usually accounted for quantitatively 



with the help of a detailed numerical parameterization (see section 3 for details). The 
second term on the r.h.s. of ( |2.22| ), i.e. 5(Ar]^ {T)M ), denotes the one-loop shift to Ar driven 
by the "genuine" (viz. non-standard) 2HDM contributions to the weak gauge boson self- 
energies. To determine this characteristic one-loop effect from the 2HDM we compute the 
set of (Higgs-mediated, h = h°, H°, A , H ± ) Feynman diagrams displayed in Fig. [j], whose 
result we indicate by Ar^-j^ , and then we subtract from it the one- loop contribution from 

H\l] 

the SM Higgs boson, H, denoted as Ar SP ^ , namely 

5(Ar£i DM ) = Ar^M " Arfif 1 , (2-24) 

see more details in the next section and in the Appendix. 

Notice that the one-loop SM Higgs contribution is included as part of the first term 



on the r.h.s. of ( |2.22| ), as we have mentioned, and it would be counted twice if it were not 
subtracted. Therefore, upon subtracting ( 2.22 ) from ( p. 16 ) and using equations ( 2.17 ) and 



(|2~2ll) and fl2~23; ), we find 



Ar 2HD M - Ar^ DM ~ Arffi? - Ar^f G) - 5(Ar^ DM ) = . (2.25) 

From now on we will identify Ar2HDM with Ar|y DM and will use Eq. ( |2.22D for the practical 
evaluation of the 2HDM contribution. 

To summarize, Eq. (2.22) encodes a reasonable estimate of the basic contributions to 



Ar in the 2HDM under the presently known perturbative unitarity constraints. It actually 
provides an upper bound to the maximum value that this parameter can reach in the 
general 2HDM after taking into account the leading one loop and higher order effects. 
This is because on the one hand it includes the corresponding SM value, and on the other 
it collects the two distinctive sources of potentially relevant 2HDM effects: i) the genuine 
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1-loop Higgs boson effects beyond the SM, i.e. <5(Ar 2 ^ DM ) ; ii) and the leading 2-loop 
0(a ew X^) effects triggered by the enhanced Higgs boson self-couplings in the limit of large 
A5, i.e. At^hdm . The latter furnishes, at large A5, the biggest source of enhancement at 
higher order within the known perturbative unitarity limits. 

Let us finish this section with a simple estimate of the numerical impact on the W- 
mass which follows from a generic shift of the parameter Ar, which we denote 5(Ar). In 
the case under consideration that shift may essentially receive the following two genuine 
2HDM contributions at one loop and beyond, which we have discussed above: 

5(Ar) = 5(ArW DM ) + Ar^ M . (2.26) 

From Eq. ( |2.3[ ) and taking the values for M% and Gf as experimental inputs our evaluation 
of Ar can be translated into a theoretical prediction for the W-boson mass, [M^y], and 
further confronted to [M^ p ]. For this we need to solve the equation 




l ™ [l + Ar(M^)] 



V2G F M% 



(2.27) 



For Ar = one obtains the tree-level value M^- ee ~ 80.94 GeV. But as mentioned above 
the full theoretical result is smaller because quantum effects imply Ar > 0. Mind that since 
Ar itself depends on M\y, Eq. ( |2,27| ) must be worked out iteratively if one aims at a precise 
prediction. To within first order, Eq. ( |2.27 ) implies that a shift S(Ar) in Ar translates into 
a small shift in the W mass given by 

1 

5M W ~--M w -2 2 6(Ar), (2.28) 
c w — s w 

with Syj ~ 0.22. In the SM, as well as in all known promising extensions of it, the quantum 
effects yield Ar of order few percent and positive, hence improving the agreement of the 
theoretical prediction with experiment. For example, in the SM Ar ~ 0.04 > 0. The 
physically measured value of M\y (M^ p ~ 80.38 GeV) is in fact smaller than the tree- level 
value by about 0.8%, i.e. some ~ 600 MeV smaller. This is consistent with (|2.28| ) since 
Ar itself is small and so that equation applies reasonably well to the entire Ar value. 

Furthermore, Eq. ( [2.28 ) tells us that even a tiny departure |5(Ar)| ~ 10~ 3 from physics 



beyond the SM induces a shift of |<5M\y| — 16 MeV in the W mass. This shift is significant 
since it is very close to the present experimental error in the W mass (M^ p = 80.385 ± 
0.015 GeV). It follows that further improvement (lessening) of that error will enable us 
to discriminate very subtle quantum effects, perhaps digging already into physics beyond 
the SM. Interestingly enough we have indicated above that the 2HDM can induce genuine 
shifts of this order from the two sources indicated on the r.h.s. of fl2.26| ). In the next 
sections we shall confirm by explicit numerical analysis that the maximum size of the total 
shift 5(Ar) can be of order of a few times 10 -3 and displace the value of M]y a few tens of 
MeV - it can reach as high as 5M\y ~ 35 — 40 MeV in the optimal cases presented here. 
This potentially significant thrust from the genuine quantum effects in the 2HDM sector 
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can bring the theoretical prediction on Mw larger and hence offers a better agreement 
with the experiment than the SM prediction, which persistently tends to stay too low as 
compared to the experimental value. 



3. 2HDM prediction for Ar and Mw at full one loop. Detailed analysis 

Hereafter we carry out a dedicated numerical analysis of the different EW precision quan- 
tities under scope. To start with, we focus on the pure one-loop evaluation of Ar, and 
thereby of M^, in the framework of the general 2HDM. Let us emphasize from the begin- 
ning the meaning of our denomination "full one loop" in the title of this section. It refers 
to the order of perturbation theory at which we compute the genuine 2HDM effects on the 
complete Ar quantity, and not just to bp. However, we incorporate in our calculation of 



Ar (following the procedure explained in Sec. 2A) the known higher order effects within 
the SM because otherwise it would be meaningless to compare with the experiment in view 
of the current high precision of the measurements. It is only in a second stage (cf. Sec. 
H]), where we perform a closer look to the genuine higher order effects emerging from the 
2HDM. Specifically we will focus there on the parameter 5p within the 2HDM and examine 
the Higgs boson self-interactions and their enhancement capabilities therein. In doing this 
we will find appropriate to resort to an effective Lagrangian approach for a first reliable 
approximation to these higher order genuine 2HDM contributions. 

The unrenormalized gauge boson self-energies Sy 101 ^ 1 , as defined in Eq. (j2.6|), reunite 
all the information we need to compute the 2HDM contributions to Ar at one-loop: in our 
notational setup (cf. Eqs. (gjp and ( Hi) ) the ful1 

one-loop payoff is denoted 



^ r 2HDM — ^ r SM + ^(^ r 2HDM) ' fi.l) 



where the second term on the r.h.s. of this equation was defined in (2.24). 

The various contributions are described by the collection of Feynman diagrams we dis- 
play in Fig. |]. For a consistent calculation of Ar 2 ^ DM we must remove the existing overlap 
between the one-loop 2HDM effects carried by the light, neutral CP-even Higgs boson (h°), 
and those driven by the SM Higgs boson (H SM ) - which are contained in ArsM itself. This 
is achieved, in practice, by subtracting from S\y,z the h°-mediated one-loop diagrams in 
the SM-like limit /3 = a — tt/2. What remains is the genuine 2HDM effect encoded in 
<5(Ar2^ DM ). Explicit details are given in the Appendix. We perform our calculation in 
the 't Hooft-Feynman gauge and regularize the UV divergences using standard dimensional 
regularization [93], which is both Lorentz and gauge invariant. Analytical manipulations at 
this stage are carried out with the help of the packages FeynArts and FormCalc [89,94]. 

Let us now describe our numerical setup. As for the Higgs boson spectra, we define 
several mass choices (cf. Table |l]) that span over representative parameter space regions, 
featuring in particular the ~125 GeV Higgs-like resonance spotlighted by ATLAS and CMS 
and in agreement with all phenomenological limits discussed in Sec. |2.2| . Set I describes a 
light, neutral CP-even Higgs state [h°] taking on the role of the 5a Higgs-candidate. It is 
thus preferably studied in the decoupling regime a = (3 — 71/2. Notice that this scenario 
can in principle accommodate a fermiophobic heavy Higgs [H°] (a = 0) for type-I 2HDM's. 



- 19 - 





M h o [GeVJ 


M H o [GeVJ 


M A o [GeVJ 


M H ± [GeVJ 


Set I 


125 


135 


200 


220 


Set II 


115 


125 


200 


220 


Ci j TTT 

Set III 


125 


200 


200 


215 


bet IV 


1 OCT 

125 


300 


300 


310 


bet V 


1 OCT 

125 


12b 


260 


300 


bet VI 


125 


126 


210 


250 


^pf VTT 


1 9^ 


1 9fi 


1UU 


9on 


Set VIII 


125 


126 


110 


150 


Set V (M H ± < M A o) 


125 


126 


300 


260 


Set VI' (M H ± < M A o) 


125 


126 


250 


210 


Set VII' (M H ± < M A o) 


125 


126 


200 


160 


Set VIII' (M H ± < M A o) 


125 


126 


150 


110 



Table 1: Higgs boson mass sets employed throughout our computation. Owing to the phenomeno- 



logical restrictions from i?-physics observables (cf. Sec. 2.2), Sets IV and V can be realized both 
for type I and type II 2HDM's, while the other scenarios would mostly be suitable for a type-I 
2HDM setup only. Sets III and IV have been specially devised to reproduce the characteristic mass 
splittings of the MSSM Higgs bosons (see the text for details). Sets V-VIII feature two neutral, 
CP-even states with masses around 125 GeV, and a constant splitting [|A/ao — A/ H ±| = 40 GeV] 
between their CP-odd [A ] and charged [H^] companions. The two possible mass hierarchies, which 
we refer to as direct [Af H ± > M^o] and inverted [Af H ± < -Mao] (sets marked with a prime) are 
analysed separately. 



Further on we shall entertain this possibility, even though a more accurate determination 
of the fermionic decay modes of the putative Higgs boson candidate might soon rule out 
this region of the parameter space. In Set II we still assume a type-I 2HDM, but in contrast 
to the previous set we identify the ~125 GeV resonance with the heavier neutral CP-even 
Higgs field [H°], whereas the lighter h° is assumed to have escaped detection so far. In this 
latter case we are of course assuming a — /3, since this choice insures SM-like couplings 
of H° to gauge bosons, and suppressed couplings for h°. Sets I and II, thus, represent 
alternative choices of Higgs bosons, featuring a non-constrained mass spectrum and in full 
compliance with the current experimental picture. Sets III and IV, on the other hand, 
assume (as in the case of Set I) a = [3 — ir/2 and are devised to mimic typical MSSM-like 
spectra - say a relatively light h° state playing the role of the ~125 GeV resonance along 
with heavier (almost mass-degenerate) companions [H^A ,^]. The mass spectrum in 
both sets has been numerically derived for selected MSSM parameter configurations with 
the help of FeynHiggs [95]. In particular, Set III corresponds to the so-called decoupling 
regime of the MSSM for maximal mixing (cf. the recent analysis of Ref. [96]). Similarly, 
we define Set IV, now with heavier masses for the scalar companions of the light Higgs 
boson, as a realization of the decoupling regime for typical mixing [96] of the MSSM. Let 
us note that owing to the phenomenological restrictions described in Sec. 2^2, Set III must 
correspond to type-I while Set IV can be either type-I or type-II 2HDM. 
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Finally, we construct the two groups V-VIII and V'-VIII' of mass sets. In both groups 
we have two neutral, CP-even, almost mass-degenerate Higgs bosons with ~ 125 GeV 
while the charged Higgs and the neutral CP-odd states are typically heavier and with a 
constant mass splitting |M^o — M H ±| = 40 GeV. The distinction of the two groups lies 
in the two possible mass hierarchies: M H ± > M^o (which we conventionally dub direct 
hierarchy); and M H ± < M^p {inverted hierarchy - corresponding to the sets V'-VIII' 
marked with a prime). The fact that the two CP-even states lie around ~ 125 GeV is 
particularly interesting, as they enable the two complementary regions a ~ (3 — it/ 2 and 
a ~ (3 to be phenomenologically viable - alongside with part of the intermediate regimes. 
This is so because, regardless of whether the trigonometric couplings are arranged as in the 
decoupling (SM-like) limit [sin(/3 — a) ~ 1]; or, conversely, if they are anchored such that 
[cos(/3 — a) ~ 1], in either case we find a neutral, CP-even state with SM-like couplings to 
the gauge bosons. The latter can thus take on the role of the ~ 125 resonance unveiled at 
the LHC - while its mass-degenerate companion will remain virtually decoupled. 



As reflected in Eq. (2.22), the 2HDM prediction for Ar is generated in our approach 
from the combination of the SM result ArsM and the genuinely new pieces sourced by the 
2HDM degrees of freedom. Since we are interested here to evaluate the full one-loop 2HDM 
correction beyond the conventional SM effects [Arl>^ DM ], we will include in this section the 
5(Ar2^ DM ) piece from ( p. 22 ) only. Recall that the other piece in ( |2.22| ), Ar^ 5 ^, first 
enters at two loops. Its quantitative impact will be assessed in the next section, as it 
requires a careful theoretical digression before we can tackle it. 

In practice we shall evaluate ArgM from Eq. (|2.3j) upon extracting the SM theoretical 



(3.2) 



value [M|^ ] from the compact numerical parametrization of Ref. [55,83,85], 



SM 

w 



Myy - di dH - d 2 dH 2 + d 3 dH 4 + d A (dh - 1) - d 5 da + d 6 dt 



d 7 dt' 



. dH dt + dg dh dt — d\Q da s + du dZ, 



The coefficients introduced thereof read 

da 



K 

di 

do 

d 3 

d,,{ 



80.3800 GeV; 

0.05253 GeV; 
0.010345 GeV; 

0.001021 GeV, 
-0.0000070 GeV, d w 

1.077 GeV, d n 



d7 
d s 
dg 



0.5270 GeV; dH 

0.0698 GeV dt 
0.004055 GeV, da 

0.000110 GeV, da s 



0.0716 GeV, 
115.0 GeV, 



dZ 
dh 



m 1 100 GeV ) 



, 174.3 GeV 
Aq 



1: 



0.05907 
a s (M z ) _ _ 

0.119 ' 



M z /(91. 1875 GeV) - 1 



100 GeV 



(3.3) 



with the top-quark mass [mt = 173.5 GeV] and the Z-boson mass [Mz = 91.1875 GeV] 
evaluated at their current best averaged values [39], while a s (Mz) = 0.118 and Aa = 
0.05911. Notice that = 80.3800 GeV in the above parameterization is not meant to 
be the zeroth order value of M\y (which we already indicated previously as M^ ee ), but a 
fiducial mass value which is chosen very close to the central value of the experimentally 
measured Jl% — confer Eq. ( |3.5| ) below. Therefore all the correcting terms on the r.h.s. 
of (|3.2D are expected to be small for any reasonable model, and some of them are very 
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much close to zero because the involved physical quantity is known with high accuracy 
(e.g. dZ ~ 0). 

The parametrization (|3.2| )-( |3~lj| ) reproduces the full set of available quantum corrections 
within the SM. It relies on the following decomposition of Ar in all the one-loop and 
higher order effects computed to date [55,83]: 

Arq M = Ar (a) + IaA^ + Ari^ + Ar (aQs) + Ar (aa ^ + Ar (G f ,a ' ra ' ) + A r ( G F m t) 

Divi i 5 os term 1 

(3.4) 

For more clarity we keep in ( |3,4j ) the notation of the aforesaid references, and hence the 
first term Ar^ denotes the complete one-loop contribution; in our notation, it corresponds 
to Ar™ in flOP . The remaining terms are higher order contributions which have been 



SM 

explicitly computed in the literature over the years 6 ; thus Ar^J [86,97] and Ar|™2 [83,97, 
98] stand for the respective bosonic and fermionic electroweak two- loop corrections; AA aaa ^ 
and AA aa ^ for the corresponding two-loop and three- loop QCD corrections [90]; and finally 
we have the pure electroweak 0(Gp mf) and mixed electroweak-QCD 0(G^a s mf) terms, 
which track the leading 3- loop contributions [99]. The remaining theoretical uncertainty 
from the unknown higher order corrections is estimated to lie around A ~ 4 MeV [38] . 
In accordance with the prescription ( p. 22 ), we then add up to the previous result the genuine 



2HDM piece 5(Ar^ DM ) from the 2HDM Higgs boson- mediated self energies [£\y,z]) as 
displayed in Fig. [l] (see the Appendix for explicit analytical details). Finally, we rephrase 



the overall Ar yield obtained in this way in terms of M$ by iteratively solving Eq. ( 2.271 ). 



This is our estimate of Ar2HDM at this point. As explained in the beginning of this section, 
we call it the "full one-loop 2HDM result", in the sense that all of the genuine one-loop 
effects from the 2HDM degrees of freedom have been taken into account on the Ar quantity 
within the 2HDM, but at the same time we have incorporated in this result the complete 
current knowledge of the leading SM quantum effects at various orders of perturbation 
theory following the above method. 

The upshot of our numerical analysis is summarized in Figs. ||-[5|. Figures |3| and [I| 
display the full one- loop 2HDM predictions for Ar and iWw, respectively, as a function of 
the neutral, CP-odd Higgs boson mass [M A o]. We sweep over a range 190 < M A o < 350 
[GeV] while varying the masses of the remaining Higgs fields accordingly, so that the mass 
splittings among the heavier 2HDM Higgs fields are mantained. At the same time we keep 
fixed the mass of h° at the value M h o = 125 GeV. For Set II we actually anchor both 
M h o and M H o at the constant values indicated in Table [l] - recall that, in this scenario, 
the heavy CP-even neutral field [H°] is meant to describe the putative ~ 125 GeV Higgs 
boson. We carry out the calculation at a fixed value of tan (3 = 1.5, in compliance with the 
stringent conditions dictated by B® — B® and B{b — > sj) physics, including the perturbative 
unitarity and vacuum stability bounds (see Sec. p.2|). On the same grounds we enforce 



6 Recall that Ar in the defining equation ( |2.3[ ) collects the expanded form of the various contributions 
up to the order of perturbation theory where the calculation has been carried out. In contrast to the form 



sketched in (2.23), no resummation is assumed here of the Aa and Sp effects since the terms obtained e.g. 
at two loop order from the resummation of leading one-loop effects are already contained in the explicit 
two-loop contribution. 
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Figure 3: Full one-loop evaluation of the parameter Ar in the 2HDM, cf. Eq. ( |3.l| ) (solid line, in 
black), as a function of the neutral, CP-odd Higgs boson mass [M^o]. The masses of the remaining 
heavy Higgs bosons (M H o, M H ±) are also varied alongside, preserving the mass splitting that defines 
each of the mass spectra (cf. Table [l]). The vertical line signals the point in the Higgs boson mass 
range that corresponds to the exact mass values quoted in Table |l|. In the particular case of Set 
II, in which the SM-like Higgs boson is identified with the heavy CP-even Higgs boson [H°], its 
mass Mjjp is kept constant, so that in practice we scan over Mpp and M H ± only. The calculation 
is performed at fixed tan/3 = 1.5, with a — f3 — ir/2 (resp. a = j3) for Sets I, III and IV (resp. Set 
II), in agreement with the SM-like decay patterns of the 125 GeV Higgs boson candidate. As a 
reference, we also display the SM prediction ArgM (dashed line, in blue) and the experimental value 
Ar exp (dotted-dashed line, in red). The associated la and 2a C.L. bands are explicitly indicated. 



SM-like Higgs/gauge boson couplings for the corresponding ~ 125 GeV Higgs boson in 
our numerical analysis - which amounts to work in the decoupling limit a = (3 — tt/2 for 
Sets I, III and IV, and a = /3 for Set II. Finally, Sets V - VIII and V - VIII', with two 
neutral, CP-even, ~ 125 GeV states, interpolate between both scenarios, and are examined 
in Figure |5[ We shall comment on them later on. 

The different panels of Fig. ||| present the Ar2^ DM results for Sets I - IV of Higgs bosons 
masses. Therewith we also display the SM value ArsM> computed as described above, 
together with the experimental value Ar exp obtained from Eq. (|2.3| ) using the experimental 
inputs 
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Figure 4: Full one-loop prediction for the W-boson mass [M$ = Af ™ ± S M^ DM ] in the general 
2HDM as a function of the neutral, CP-odd Higgs boson mass [Mao]. The results we present as a 
deviation [SM W ] = M$-M^ p with respect to the experimental value M^ p = 80.385 ±0.015 GeV. 
The masses of the remaining heavy Higgs bosons (M H o,M H ±) are also varied alongside, preserving 
the mass splitting that defines each of the mass spectra (cf. Table |l|). The vertical line signals 
the point in the Higgs boson mass range that corresponds to the precise mass values quoted in 
Table [l]. In the particular case of Set II, in which the SM-like Higgs boson is identified with the 
heavy CP-even Higgs boson [H°], its mass Mho is kept constant, so that in practice we scan over 
Mao and M H ± only. The calculation is performed at fixed tan/? = 1.5, with a = f3 — it/2 (resp. 
a = f3) for Sets I, III and IV (resp. Set II), in agreement with the SM-like decay patterns of the 125 
GeV Higgs boson candidate. The green band provides an estimate on the theoretical uncertainty. 
As a reference, we also display the SM prediction M|y (dashed line, in blue). The associated la 
and 2cr C.L. bands for the measured value M^ p are explicitly indicated. 



M exp 
1V1 W 



80.385 ± 0.015 GeV M z -- 
a(0) = 1/137.03599 G F 



91.1876 ±0.0021 GeV 
= 1.16637 10- 5 GeV" 2 . 



(3.5) 



Plugging the above values into Eq. (2.3) we get 



Ar, 
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V2G F 
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w 



1 = 36.322 x 10" 3 . 



(3.6) 
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Figure 5: Full one-loop evaluation of the parameter Ar in the general 2HDM, cf. Eq. (3d) (solid 
line, in black) as a function of sm((3 — a). The calculation is performed at fixed tan/3 = 1.5, for 
sets V to VIII and V to VIII' of Higgs boson masses, and for the two different Mao /M h ± mass 
hierarchies. As a reference, we also display the SM prediction ArsM (dashed line, in blue) and the 
experimental value Ar cxp (dotted-dashed line, in red). The associated la and 2a C.L. bands are 
explicitly indicated. 



In Fig. |] the corresponding 2HDM and SM theory predictions for M$ are represented as 
a function of M^o under the same conditions as in Fig. || One can also appreciate the 
deviations of the theoretical W-boson mass with respect to the latest experimental value 
M^ p = 80.385 ± 0.015 GeV. Confidence levels for M^ p and A 

'"exp are also included in 

our plots as colored areas corresponding to the la (yellow) and 2a (magenta) regions. 
In the case of Ar, we compute these uncertainties from those of by standard error 

propagation methods. The region comprised between the two parallel thin strips (in green) 
in Fig. [I] supplies an estimate of the theoretical uncertainty [A M^]. The latter we quantify 

by means of Eq. ( 2.27 ) upon substituting l + 5(Ar^ i - DM ) — > 1/(1 — ^A^hdm))' wmcn leads 
to a slightly modified prediction [M] . Such a shift operates an approximate resummation 
of the 2HDM contributions, in such a way that the difference relative to the central value, 
[AM$/M$ = \M$ - M$|/M$], should be indicative of the size of the higher-order 
effects neglected when truncating our perturbative expansion at 0(a ew ). 

Complementary vistas on the Ar behavior across the 2HDM parameter space are 
displayed in Fig. ||. We compute once again the full-fledged one-loop 2HDM prediction 
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[Ar2y DM ] and compare it to the SM result [Ar$M] and the experimental value [Ar cxp ]. 
We keep a fixed tan/? = 1.5 but we now sweep all over the sin(/3 — a) range. In doing 
so, we effectively interpolate between the two corners with SM-like Higgs/gauge boson 
couplings. In the light of the LHC findings, regions away from sin(/3 — a) ~ 1 - or 
conversely sin(/3 — a) ~ - would be disfavored, if not simply excluded, at a certain 
confidence level. We will refrain from considering a more accurate treatment of this issue, 
as it is still subdued by large statistical uncertainties and surely not relevant for our present 
discussion. Fig. ||| displays a featureless Ar2HDM profile as a function of sin(/3 — a), which 
follows from the approximate mass-degeneracy of the neutral CP-even Higgs boson states. 
Worth noticing is also the influence of 5p (what is tantamount to say, the size of Ar^ p l in 
the notation of Eq. ( 2,11| )): the lighter the Higgs boson masses -and the narrower their 



splittings - the smaller becomes Ar^. One can easily understand this behavior from 



equations ( 2.15 ) and ( |6.6[ ). One visible signature that stands out in these plots is that the 
tension between the theoretical prediction [A^hdm] and the experimental value [Ar exp ] 
shrinks significantly when comparing the outcomes for the different mass choices from 
Sets V to VIII. For Set V, for instance, we have \5p\ ~ 0(1O -3 ) - at the very border of 
the permitted custodial symmetry bounds. The predicted Ar2HDM is then pushed more 
than 2a away from the experimental result - with an even larger tension than the SM 
prediction. In contrast, for Set VIII the theoretical prediction [A^hdm] moves closer to 
[Ar exp ], departing from it roughly ~ 0.6% - the SM prediction staying circa ~ 4% away. 
All these effects are qualitatively similar, although opposite in sign, when we swap the 
mass hierarchy. Remarkably, we find that [Mao > Af H ±] pulls our prediction for M$ 
slightly further away from M^ p , as compared to [M^p < M H ±]. For the scenarios we have 
considered, therefore, we can conclude that the EW observables we consider here tend to 
favor 2HDM realizations with the charged Higgs being heavier than the neutral, CP-odd 
one. 

The predicted M$ values within the 2HDM are seen to follow a smooth and monotonous 
variation, with changes of O(10) MeV - equivalently, 5(Ar) ~ 0(1O -3 ) as retrieved by the 
previous Fig. || - when sweeping the Higgs boson mass range 190 < M A o < 350 [GeV]. 
The largest possible theoretical W-boson mass shift from genuine 2HDM effects is sub- 
stantial: \5Mffi DM \ ~ 35 MeV ( more than twice the current experimental error on M\y). 
Larger variations ensue from a trade-off between two different - and somewhat opposite - 
conditions, which are the following: i) heavier Higgs bosons, in correspondence with the 
mass suppression of the one-loop Higgs- mediated contributions to Sw,Z) yield relatively 
tamed quantum effects; and ii) broader splittings among the different Higgs boson masses, 
lead to a stronger breaking of the approximate SU(2) custodial symmetry, and hence to 
an augmented contribution Ar^. A similar balance between both tendencies we will also 
identify when examining the leading higher order quantum corrections, as discussed further 
down. 

The SM theoretical prediction following from the full parametrization ( |3.2| )-( |3l3| ) with 
M H = 125 GeV reads 

M^ M = 80.363 GeV. (3.7) 
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The difference M^ — M^ p corresponds to the horizontal line marked as SM in all the plots 
of Fig.||] - except for Set II, for which the identification h° = H implies that M\y is to be 
evaluated at M H = 115 GeV, yielding M w = 80.367 GeV. We see in these plots that the SM 
value is systematically too low, it typically fails to match the central value Af^ p by more 
than 20 MeV below it. This deviation is not worrisome right now because it lies within the 
la region, but this could change in the future when the error in M^ p may dramatically 
decrease at the ~ 5 MeV level. At this point the SM prediction, whose intrinsic theoretical 
error is estimated to be ~ 4 MeV [38], could be in good shape or on the contrary it might 
already clash with experiment. In contrast, the 2HDM numerical prediction is seen to have 
the ability to easily cure the differences with respect to the experimental value, should the 
conflict with the SM prediction become confirmed. From Fig. Q it is evident that the 2HDM 
prediction can stay very close to the central experimental value M^ p and moreover it can 
accommodate deviations of up to 20 — 35 MeV with respect to it, depending on the Higgs 
mass sets used from Table |l[ 

Let us remind the reader at this point that within the approach employed in this 
section, we are still not sensitive to the potential effects from the A5 parameter. Neither 
the results show any direct sensitivity to whether they are obtained within a type-I or a 
type-II 2HDM setup, as this depends on the specific Yukawa couplings of each model. Still, 
we should bear in mind that both features, namely the strength of the Higgs self couplings 
and the pattern of Higgs /fermion Yukawa interactions are indirectly present through their 
interplay with the theoretical and phenomenological constraints - which are duly taken 
into account in our analysis (cf. Sec.|2.2|). Let us recall in particular the lower bound of 
~ 300 GeV for the charged Higgs boson mass for type II models. 

We conclude this section by emphasizing what we believe is our most important obser- 
vation, hitherto unnoticed (to the best of our knowledge) in the literature, to wit: the bulk 
presence of the term <5(Ar 2 y DM ) - defined in Eq. ( 2.24[ ) - provides a significant shift on the 



VF-mass genuinely caused by the 2HDM heavy Higgs companions. This shift softens the 
influence of the SM(-like) Higgs boson [h°] on the W-boson mass prediction and therefore 
relaxes the current M^-versus-M^ p tension, which is well known to grow with M h o and 
constitutes one of the main reasons why global fits to Electroweak precision data tend to 
favor light Higgs bosons - rendering the minimum x 2 -values for masses even below the 
LEP bounds (cf. e.g. [101]). What our results reflect is that such trend no longer holds in 
the context of the general 2HDM. Here, and essentially for all the surveyed scenarios, the 
neat effect of adding new scalar SUl(2) doublets leads to a systematic downward shift of 
Ar, and so upward on M$, which tends to improve the agreement with the corresponding 
experimental measurement. Said differently, when compared to the SM expectations for a 
~ 125 GeV SMdike Higgs boson mass, the genuine 2HDM contribution [(5(Ar 2 ^ DM )] tends 
to pull down the overall Ar 2 y DM prediction. Correspondingly, M$ becomes larger, and 
hence closer to the corresponding experimental measurement M^ p . We deem that this 
is the most relevant message to be conveyed at this point of our analysis. In the next 
section we focus on the additional influence of the 2HDM quantum effects through the A5 
parameter. 
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Figure 6: Pictorial sketch of the effective Lagrangian approach devised to capture the leading 
2HDM higher order corrections to 5p. In the limit of large Higgs boson self-interactions, these 
are dominantly driven by the Higgs- mediated contributions to the weak gauge boson self energies, 
as displayed in the lower part of the figure (cf. Fig.|| for a specific sample). As a first stage we 
compute the one- loop Higgs- mediated corrections to the different Higgs/gauge (and Goldstone) 
boson interactions to order (9(Af); these correspond to a UV finite, gauge invariant subset of 
contributions which encapsulates the (dominant) effects triggered by the enhanced Higgs boson 
self-interactions (upper- left corner). Such corrections we then rewrite as one- loop form factors and 
use them to effectively dress the leading-order Higgs/gauge boson couplings (upper- right corner). 
The resulting set of Born-improved effective vertices (depicted as colored blobs) can be finally 
plugged back in the calculation of the one-loop gauge boson self energies (lower-right and lower-left 
corners). 



4. Higher order 2HDM effects from enhanced Higgs boson self-interactions 

A number of recent studies [64, 81] have pointed out and exploited the possibility that the 
Higgs/gauge boson couplings eventually undergo a drastic modification at the quantum 
level in the limit of large Higgs boson self-interactions. These modifications are in practice 
restrained by certain theoretical and phenomenological bounds, but their net effect can 
still be quite sizeable in different contexts. In the present case the potential enhancements 
we can actually trade through the aforementioned O(X^) terms in the Higgs boson self- 
coupling structure. Such scenarios are not only phenomenologically viable, but also very 
well supported on theoretical grounds. As a matter of fact, the rise of scalar resonances 
constitutes a prominent manifestation of models with strong dynamics in the EWSB sec- 
tor [15,102]. Interestingly, strong self-interactions among Higgs, or Higgs-like states, are 
most often identified with rather heavy mass spectra - reflecting the relations among masses 
and self-couplings settled through the Higgs potential. Nevertheless, the 2HDM comprises 
a wider range of possibilities - as not all of the Higgs self-couplings are tied to the corre- 
sponding physical Higgs masses. The fact that we can freely dial the parameter A5 enables 
to recreate genuine 2HDM regimes with relatively light, and yet strongly coupled Higgs 
bosons. 



As we have preliminarily discussed in Sec. 2.4, EW precision quantities - most signif- 
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icantly Ar and 5p - become sensitive to the Higgs boson self-couplings via 2-loop effects 
on the gauge boson self-energies. This is illustrated in Fig. ||. Here we do not target 
at an explicit 2-loop calculation. Instead, we aim at getting insight into the foreseeable 
quantitative impact that such Higgs boson self-couplings may exert on the theoretical pre- 
diction for the W-boson mass [M<^]. In fact, since a fully-hedged 2-loop calculation is too 
cumbersome and probably unnecessary at this point, we resort to an effective Lagrangian 
approach. The advantage of it is that we can focus more easily on the main potential 
sources of enhancement and moreover we can even have access to their influence beyond 
the 2-loop level. We can say, in colloquial terms, that by using an effective approach we 
can better avoid not seeing the forest for the trees. 

In practice we devise our effective Lagrangian approach as follows: 7 : 

1. To start with, let us recall that after SSB of the EW gauge symmetry we can always 
trade the original Higgs self-couplings Aj (except A5) for the physical masses, mixing 
angles and gauge couplings. This is explicitly done e.g. in Sec. Ill of Ref. [64] 
- cf. Eq. (12) there. From here we can easily implement the limit of large Higgs 
self-couplings by projecting out the As-dependent part of each of the Higgs self- 
interactions in the model and neglecting the other terms. According to this procedure 
we can redehne Aj — > Aj for all the original self-couplings Aj in the Higgs potential 
-cf. Eq. Q ~ as follows: 



Ai = ^As(l-tan 2 /3); A 2 = - A 5 (1 - cot 2 /3); A~ 3 = --A 5 ; A 4 = A 6 = . 

(4.1) 

The above approximation has its counterpart in the physical basis. Following the 
splitting procedure ( p.20|) , we collect the A5 terms associated to the Higgs self- 
couplings in that basis. They render structures of the sort: 

f.6M w s w cos(a + /3)cos 2 (/3-a)\ 

A h o h o h o c h A 5 = 1 i ^— -i '- x A 5 , 4.2 

V e sin 2/J / 

and analogously for all the self-interactions Ahhh relating the different Higgs boson 
fields to each other (we refer once more e.g. to Table II of Ref. [64]). Notice that the 
strength of all the Ahhh couplings, and so also their related potential enhancements, 
are now fully determined as functions of A5 and tan (5 - with no remaining dependence 
on the Higgs nor the gauge boson masses. As it will become clear soon, this is 
a necessary step in order to preserve gauge invariance within our framework. In 
Fig. [7| we illustrate how this approximation performs when compared to the complete 
analytical form of the triple (3h) self-interations (see e.g. Table III of Ref. [64]). All 
couplings in this figure are normalized to the value of the 3H self-coupling in the SM, 



7 Approaches along these lines are certainly not foreign in the literature, see for example Ref. [27]. In 
our opinion, however, scarce attention has been devoted to the underlying assumptions and corresponding 
limitations. 
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Figure 7: Coupling strength [Ahhh] for a sample of triple Higgs self-interactions, as a function 
of the parameter A5. The couplings are normalized to the strength of the 3H self-coupling in 
the SM, assuming that the SM Higgs boson mass yields Mh = 125 GeV. As a result we get 
a hhh( a/ h = 125 GeV) = 188.24 GeV. We can check that the obtained Ahhh values lie systematically 
below the Lee-Quick-Thacker bound A|^ H (M H ~ 1 TeV) ~ 1.3xl0 4 GeV [91]. The figure compares 
the size of these 3h self- couplings when considering i) the projection of the ©(As) components only 
(left panel) and ii) its exact analytical form (right panel) - see the text for details. Without loss of 
generality, we dwell here on set I of Higgs boson masses and we fix tan/3 = 1.5. The shaded regions 
on the left (resp. right) side are excluded by unitarity (resp. vacuum stability). 

for a Higgs boson mass of M H = 125 GeV, see Eq. (|2T9|) . Noteworthy is that the 3h 
enhancements within the 2HDM, as indicated in Fig. [j] for values of | A5 1 ~ O(10), are 
perfectly consistent with the Lee-Quick-Thacker unitarity bound for the 3H coupling 
in the SM [91]: 

A _ 3eMg - 
2sw M\v - 

The main virtue of Fig. |7] is that it makes apparent how both descriptions of the Xhhh 
couplings — viz. the full one versus the truncated approach Ahhh — > defined by 
( Ug) — nicely agree in the purported limit of large Higgs boson self-interactions. 
At the end of the day, the goodness of this approximation will depend on the extent 
to which this limit is effectively realized by a given choice of A5 and tan /3. 

2. Second, we examine the leading quantum effects on the different Higgs/gauge boson 
couplings in order to absorb them into effective vertices. These we can generically 
sort out into three categories: 

i) Higgs/Higgs/gauge boson couplings: <?hhv 

ii) Higgs/gauge/gauge boson couplings: <?hvv 

iii) Higgs/gauge/ Goldstone boson couplings: <?hVG 

with [h = h ,H°,A°,H ± ] standing for a generic 2HDM Higgs field; V = [W ± ,Z°] 
for a weak gauge boson; and G = [G ± ,G°] for the associated charged and neutral 



~ 1.3 x 10 4 GeV. (4.3) 

M H =lTeV 
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Figure 8: Sample of Feynman diagrams that account for the 0{\\) one- loop corrections to the 
h°A°Z n and h°W ± H ± gauge couplings in the 2HDM. 



Goldstone modes, namely the longitudinal components of the gauge bosons in the 't 
Hooft-Feynman gauge. In practice, structures of type iii) are formally equivalent to 
those of type i), as they involve two scalar and one vector fields. In the regimes of 
interest here, namely in the large Higgs self-interaction limit, the leading corrections 
to the above couplings are thus driven by the interchange of virtual Higgs bosons. 
A sample of the corresponding Feynman diagrams is provided in Fig. ||, in which 
we single out the case of the g^pppz anc ^ fi , h°w ± H ± effective couplings. The one-loop 
corrections to each of these couplings involve both triangle diagrams and self-energy 
insertions; the former account for the genuine 0(A§) one-loop vertex corrections, 
whereas the latter involve also C(A|) pieces from the finite wave function renormal- 
ization of the external Higgs boson legs. Starting from the 1-loop vertices of Fig. ^ we 
can easily recognize them as being part of some of the vertices in the 2-loop diagrams 
of Fig. [2]. When the vertices in the former collapse to a point we obtain the kind 
of solid (colored) blobs indicated in Fig.||. These blobs represent static form factors 
attached to the corresponding vertices, i.e. evaluated at zero momentum, and are 
therefore appropriate to estimate the higher order effects on the <5/>parameter (even 
if they are actually not suitable for addressing the corresponding effects on the full 
Ar parameter, see later on). 



The practical recipe in our approach is now clear. Based solely on power counting 
and dynamical considerations, the one-loop effects from Fig.|8| can be conveniently 
cast as form factors of the guise: 
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flhhV ~ 16^2 /hhv(//M|) (p • e) ~ C(a ew A 2 ,); 

6hVV ~ 1^ (^|f) 2 ^VV(P 2 /M, 2 ) (e • e) ~ 0(a ew A 2 ); (4.4) 

ChGV ~ ^2 (^|f) 2 hGv{p 2 /M 2 h ) (p • e) ~ O(o ew A 2 ) . 

These form factors carry the sought-for fingerprint of the Higgs boson self- interactions, 
which we generically flag here as [Ahhh]- The typical Higgs boson mass scale Mh in 
Eq. fl4.4|) restores the corresponding dimensions and accounts for the correct decou- 
pling properties. At the same time, it partially balances the enhancement effect 
of the Higgs self-interactions in the numerator. By /(£> 2 / m D we denote a generic 
rational function, depending on the (ratios of the) relevant scales involved in the 
Higgs boson- mediated one-loop diagrams, namely the masses of the virtual particles 
and the transfered momentum. These form factors we evaluate at zero momentum. 
Put another way, we reduce the loop structure to a point-like interaction in which 
tfihh ~~ ^ a PP ears as the only effective coupling in the large A5 limit. Finally, 1/167T 2 
stands for the usual numerical factor from the one-loop integrals. We carry out this 
calculation with the help of FormCalc [94]. Explicit analytical details we provide 
in the Appendix, cf. Eqs. ( |6.15| )- (|6.19 ). Equivalent form factors could of course be 



entertained for the MSSM case. However, in stark contrast to the general 2HDM, 
they would show a pure 0(g 3 ) gauge structure and would therefore be inconspicu- 
ous from a phenomenological viewpoint, in the sense that no characteristic MSSM 
signature could easily emerge from them. 

Noteworthy is also the following. The one-loop form factors (4.4) computed this way 



are UV finite and gauge invariant. The finiteness follows by simple power counting. 
The gauge invariance ensues from the fact that we consistently retain all 0(a ew \1) 
contributions - but not more. It should be clear that it is because of our exclusive 
selection of the characteristic A 2 , tag in these diagrams that we are granted to suc- 
cessfully isolate a meaningful (gauge invariant) and finite contribution. By the same 
token the sum of all the original 2-loop diagrams as such carrying the A 2 , tag (cf. 
Fig.|| for a representative sample) is gauge invariant and finite. Had we kept also 
the full dependence of the trilinear couplings on the Higgs and gauge boson masses, 
we would have been forced to take into account the remaining set of contributions 
not only from the Higgs bosons, but also from the gauge sector (viz. the gauge and 
Goldstone boson-mediated diagrams) in order to insure the overall gauge invariance 
of the Higgs/gauge form factors (|4.4[). The fact that we single out these C(A|) pieces 
from the 3h self-interactions could be viewed as if we were effectively decoupling the 
dynamics of the Higgs and the gauge sectors of the 2HDM. Of course this is not a 
decoupling in the usual sense, namely it is not due to the presence of a hierarchy of 
masses - notice, indeed, that the mass spectra under examination typically feature 
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My < Mh- Rather than "decoupling" we are "detaching" the As-effects from the 
rest of the quantum contributions in a consistent way. There is actually some inde- 
pendence at the level of QFT renormalizability between, say, the gauge and Yukawa 
sectors on the one side, and the As-structure of the Higgs boson self-interactions on 
the other side. These sectors contain independent parameters (gauge and Yukawa 
couplings in the the former versus trilinear couplings in the latter). They all carry 
some characteristic renormalizability "flagpoles", which can nevertheless be mixed. 
For example, in the case of the Higgs self-interactions there is one part that still 
carries the gauge flagpole, but there is another which does not (the As-part). It is 
precisely this part that has no crosstalk with the others since it is gauge invariant 
and finite per se. As a result there can be, in principle, an arbitrary hierarchy of 
strengths between the Higgs self-couplings and the (merely gauge) Higgs/gauge and 
Higgs/fermion interactions. This is what allows to detach one sector from the other 
without jeopardizing in any possible way the QFT renormalizability of the theory. In 
turn, this is what authorizes the 2HDM Higgs bosons to couple to each other much 
stronger than they do to the gauge bosons. Such property is impossible in the MSSM, 
where the gauge flagpole is inherent (owing to the underlying supersymmetry) to the 
entire structure of the Higgs self- interactions. 

To summarize: in our calculation the A5 parameter plays the role of a flagpole which 
marks an independent sector of the theory. By selecting it we can factorize a gauge 
invariant and UV-finite piece from the entire calculation. In the limit of large A5 
such well-defined piece encapsulates the dominant quantum effects on the tree-level 
Higgs/gauge interactions. 

3. Finally, we can fold these tree-level Higgs/gauge interaction Lagrangians with the 
form factors [«hhV) ^hw> ^hcv] derived above and build up a corresponding Born- 
improved set of effective Lagrangians: 



£hw = Shhv (1 + ahhv) [(dph) V^h-hV^ d^h] ; 

£wv = <7hw (1 + khvv) h ; (4.5) 

£wg = SWG (1 + c hVG ) [(d„h) V»G-hV» d^G] . 

The general structure of the associated form factors can be illustrated e.g. in the 
case of the effective h°A°Z° coupling, 



a h o A o z o 



P 2=0 = ^ Vh ° A ° z " ~ S ' h ° ~ S A° ~ tan ^ 2p °^ Ke ^h°H°l ( 4 - 6 ) 



p 2 =0 

where the explicit form of the two and three point functions is given in the Appendix. 

The idea of absorbing the bulk of the quantum effects into improved leading-order 
effective interactions ( ft.5| ) is of course a well-known procedure within the time-honored 
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Figure 9: 2HDM contribution to the parameter Sp, as a function of the self-coupling A5, for 
different sets of Higgs boson masses. We superimpose the pure one-loop [<5/?W] and the improved 
one-loop [(5pl 1_eff l] results. The latter trades the dominant higher order effects from the enhanced 
Higgs boson self-interactions. The quantitative impact of these higher order corrections on the 
prediction of the W-boson mass is displayed in the lower panels. The latter we compute via the 
contribution from 8p^~ eS ^ to Ar, see Eq. (4.8). For each set of Higgs boson masses, the value of 
tan/3 is chosen such that the 3h self-couplings maximize for the largest allowed | A5 1 values. The 
shaded areas account for the unitarity (left) and vacuum stability bounds (right). 



effective Lagrangian approach to QFT. Here we just mention that it has been applied with 
success in the past to analogous contexts. A well known example is the enhanced bot- 
tom/sbottom Yukawa coupling in the large tan/3 limit of the MSSM [104] - see also [103]. 
The approximation will be meaningful as long as we probe these effective couplings at 
momentum scales significantly smaller than the typical Higgs boson masses. For quantities 
which are momentum-dependent, though, such effective coupling treatment may no longer 



-34- 





( T 

oet l 

Ct — p — 7T / Z 


Q-vf T 

— n 

Oc — u 


Q,4- TT 

oet it 

Ct — p 


( TT 

oet ii 

SM TT /9 

Ct — 7T / Z 


Q„j- TTT 

oet 111 

Ct — p — 7T / Z 


oet iv 

Ct — p — 7T / Z 


t-±T SM [X 1U 


*}7 7/1/1 


^7 7/1/1 


^7 AA'X 


37 AM 


37 7/1/1 


*37 7/1/1 


Ar 111 fxl0~ 3 l 


36.701 


36 650 


36 352 


36.409 


35.845 


35 603 


5(Ar^ DM )[xlO- 3 ] 


-1.043 


-1.094 


-1.091 


-1.034 


-1.899 


-2.141 


^(MiL M )/Ar^ DM | [%] 


2.8 


3.0 


3.0 


2.8 


5.3 


6.0 


<5p 111 [xlCT 4 ] 


-2.249 


-2.192 


-2.325 


-2.387 


-0.302 


-0.134 


5p ll - ctt| [xlO" 4 ] 


-2.403 


-2.348 


-2.709 


-3.014 


-0.197 


-0.160 


|A(<5p cft )/5pW| [%] 


6.8 


7.1 


16.5 


26.3 


35.0 


19.0 


5M§ [MeV] 


17.257 


18.078 


18.021 


17.102 


31.071 


34.967 


AMjJ" clt| [MeV] 


-0.865 


-0.880 


-2.166 


-3.537 


0.596 


-0.144 


tan f3 


1.36 


1.36 


1.36 


1.36 


1.38 


1.17 


A 5 


-8.68 


-8.68 


-8.68 


-8.68 


-8.46 


-11.08 



Table 2: Detailed numerical analysis of the different electroweak quantities under survey, these are: 
Ar, Sp and Mw- We examine different choices of Higgs boson masses (cf. Table |]) and trigonometric 
couplings. The values of tan/3 and A5 (see the bottom rows of the table) maximize the enhanced 
higher order effects induced by the Higgs boson self-interactions. The notation S(Ar^ BM ) spells 
out the "genuine" one-loop 2HDM effects (i.e. after consistent subtraction of the SM part) - cf. 
Eqs. fl2.22| ) - fl2.24| ). In turn, 5p^ and 5 M$ denote the one-loop shifts on these parameters from 

those genuine 2HDM one-loop effects, while A(Sp cS ) and AM£~ ei ' represent the corresponding 
higher order effects on Sp and Mw beyond one-loop, computed in our effective Lagrangian approach 



cf. Eqs (4.5) and (4. 



be suitable. But in our case the approach is fully justified, as we are using this approxi- 
mation only to estimate some higher order effects on the 5p parameter, a static prominent 
piece of Ar, namely Ar^ = — (c 2 ^/s 2 v )5p, for which the gauge boson self energies are 
to be evaluated at zero momentum - cf. Eq. (^^). This is the reason why the effective 
coupling treatment elaborated here is not directly exportable to e.g. a full-fledged calcu- 
lation of Ar. As we will discuss in more detail, the Ar^' component proves insufficient 
to completely describe the quantum effects from the Higgs bosons, which means that the 
momentum-dependent part of the gauge boson self energies (evaluated on-shell p 2 = My) 
cannot be neglected, and in fact may be quantitatively relevant. 

With these effective Lagrangians at hand, an estimate of 5p beyond one-loop - viz. 
including the dominant 0(a ew A|) terms - is now at reach with a moderate, albeit still non- 
trivial, amount of work as compared to the full 2-loop calculation. The obtained numerical 
results are presented in the upper panels of Fig. ||. Here we plot the evolution of 5p versus 
the Higgs self-coupling A5 and for the Sets I- IV of Higgs boson masses quoted in Table [j]. 
For each set, we fix tan /3 so that the quantum effects governed by the enhanced 3h self- 
interactions maximize for the largest attainable |As| values within bounds. Alongside the 
pure 2HDM one-loop effects [<5/?M = dp^^y,] - which correspond to the flat, dashed line 
- we superimpose the improved, A5-dependent, curves 

[fy[l-eff]] derived via the effective 
Lagrangian interactions defined by Eqs. (|4.5| ). By dialing the parameter A5 we portray the 
dependence of 5p^~ eS ^ on the Higgs boson self-interaction enhancements. With decreasing 
A5 values, the difference A(5p eff ) = (5p^~ eS ^ — 5p^) obviously goes to zero - and in this 
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Table 3: Coupling strengths Ohhv an d &hvv ~ cf. Eq. (4.4) - for a representative set of hhV and hVV 
effective interactions. We compare them numerically for Sets V and VIII (cf. Table |]), assuming 
a = /3 — it/2 and fixing tan/3 and A5 such that the triple Higgs boson self-coupling enhancements 
maximize - in a way compatible with all the bounds. 



case the curves bp^ and bp^~ eS ^ (A5) tend to merge. 

A complementary numerical account is provided in Table ^. Here we consider Sets I- IV 
of Higgs boson masses from Table p], with corresponding SM-like couplings for the ~ 125 
GeV CP-even Higgs boson or H°) state in all the cases. Moreover, we also investigate 
the particular instances a = (for Set I) and a = ir/2 (for Set II), these are the so-called 
fermiophobic Higgs limits of type-I 2HDM. Most significantly, in this table we quantify 
the maximum attainable departures |A(£p )| from the plain one-loop predictions. These 
optimal scenarios are realized for particular choices of tan/3 and A5 (as quoted in the last 
rows of the table), saturating the 2HDM unitarity bounds [92]. Such configurations are 
preferably achieved for small tan/3 ~ 1 — 2 and moderate (negative) A5 ~ —5/ — 10. 

These results clearly spell out the physical meaning of bp as a measure of the SU(2) 
custodial symmetry violation. The plain one-loop prediction for bp ranges from O(10~ 4 ) 
(for Sets I and II, viz. for relatively light Higgs bosons with unconstrained mass splittings) 
to 0(1O -5 ) (for Sets III and IV, wherein the heavier Higgs bosons are tailored to mimic 
the mass splittings of SUSY-like spectra). This same feature we may track down too from 
the right panels of Fig. |l(], in which we again superimpose the mere one-loop prediction 
[JpM] and the improved one [Jp[ 1_eff l], as a function of A5, for the complementary mass 
Sets V-VIII. The sizable value of bp, at the 0(1O -3 ) level for Set V (this means nearly 
overshooting the custodial symmetry limit \bp\ < 10 -3 ) depletes very remarkably when the 
heavier Higgs boson masses - and so the different mass splittings - are pulled down all the 
way from Set V to Set VIII (cf. Table §. 

By the same token, the relative importance of the leading higher order corrections to 
bp becomes drastically promoted in some cases (at the highest available values of A5), as 
we quantify explicitly on the left panels of Fig. |l0[ This is in part because these higher 
order quantum effects are added to a decreasing one- loop piece [bp^] - which falls down 
by roughly one order of magnitude if we compare again Sets V and VIII. At the same 
time, the boost in the relative departure A(bp eS )/bp^ is also partly explained due to 
the smaller mass suppression of the one-loop Higgs-mediated effects. In Table || we settle 
this statement quantitatively, by comparing the actual coupling strength of a number of 
effective interactions for Sets V and VIII - this is to say, when we assume heavier (resp. 
lighter) Higgs bosons (cf. Table [j]). 
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Figure 10: Impact of the leading higher order corrections to Sp. On the left panels we evaluate their 
influence through the relative departure [A(<5 efI )/<5pW = (Sp^^ cti ^ — Sp^)/Sp^ [%]] as a function 
of the self-coupling A5 . We consider the Higgs boson mass sets from Table [l] with both direct (Sets 
V-VIIL top panel) and inverted (Sets V'-VIII', bottom panel) mass hierarchies between the A and 
the fields. We set a — (3 — ir/2 and fix tan/? as quoted in Table [[[ so that the 3h self-coupling 
enhancements maximize for a given A5 - saturating the unitarity bounds for each of the mass sets, 
cf. Table ^. Complementarily, on the right panels we present the absolute value of 5p, again as 
a function of A5 and for the different sets of masses. In each case we superimpose the mere one- 
loop results [<5pM] and those including the dominant higher-order corrections [<5p[ 1— cff J (A5)] . Let us 
notice that unitarity and vacuum stability constraints are not explicitly included in the figure; the 
phenomenologically viable A5 range for each of masses is explicitly quoted in Table ||. 



We can also rephrase these results as a relative contribution to Ar (cf. Eq. ( 2.11 )), 
namely via the term Ar^l, of the order of ~ 5% to ~ 50% to the overall 2HDM prediction 
[A^hdm]- This quantity grows along with the mass splittings between the different Higgs 
bosons. That explains e.g. why, for Sets I and II, these <5p-driven contributions are more 
sizable than those for the SUSY-like Sets III and IV (that stagnate at the 0(1O -5 ) level). 
In contrast, for the last two sets the total Ar values are larger. This simply means that the 
on-shell gauge boson self-energies or, in other words, the non-static contributions to Ar, 
are more relevant here, as compared to the self-energies evaluated at zero momentum. The 
induced shift (from strict one- loop 2HDM effects) on the W-boson mass prediction [AM^] 

is correspondingly larger for the SUSY-inspired mass sets. Specifically, AM$ ~ 15MeV 
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Figure 11: Overall shift in the theoretical prediction for the W-boson mass (cf. Eq. ([Of)), 
driven by the leading Higgs-mediated higher-order corrections, as a function of the self-coupling 
A5. The results we display separately for sets V-VIII (cf. Table [l]) with a direct mass hierarchy 
[M H ± — M A o > 0] , and for different choices of the mixing angle a. We fix tan (3 as quoted in Table ^, 
so that the 3H self-coupling enhancements maximize at the largest allowed |A 5 | values - also given 
in the table. 



for Sets I and II, versus AAf^ ~ 30 MeV for Sets III and IV - cf. Table | and Fig. £§ This 
is not surprising, after all. It is well-known that the MSSM tends to relax -sometimes very 
substantially - the existing tension between the SM theory predictions and the experimental 
measurements of electroweak precision quantities - for a recent analysis see e.g. [105]. We 
may turn around this observation and say - emphasizing once more the claim made in Sec. 
I - that our study shows that the general 2HDM can perform equally well than the MSSM 
in the task of fostering the agreement between the theoretical prediction of the W-vaass 
versus its experimentally measured value. 

In this respect, it is interesting to point out that in the MSSM case the Higgs bosons 
alone induce - if barring a very heavy, and hence virtually decoupled SUSY spectrum - 
a rather moderate contribution to the EW quantities Ar and 5p, when compared to the 
squark, slepton and chargino-neutralino sectors. Let us be more quantitative here. Set III 
of Higgs boson masses (cf. Table [l]) is built upon a specific MSSM configuration whose 
corresponding SUSY spectrum consists of 0(100) GeV charginos and neutralinos, together 
with ~ 700 GeV stops, large Higgs-stop trilinear couplings, and the remaining sfermion 
masses lying above 1 TeV. Employing FeynHiggs to extract the full MSSM contribution 
to 5p within this scenario, we find <5/9mssm = 2.186 x 10 -4 . In addition we can separately 
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evaluate the contribution driven by the MSSM Higgs sector alone by means of our 2HDM 
calculation, namely by plugging the mixing angle value «mssm = —0.188 which corresponds 
to that specific MSSM benchmark. Once more we can use FeynHiggs for the numerical 
evaluation, and we obtain ^/O^sIm = —2.910 x 10 -5 . This result clearly illustrates how the 
pure SUSY Higgs-induced contribution to 5p falls roughly one order of magnitude below 
the total MSSM budget. The SUSY embedding for Set IV, in turn, features much heavier 
[viz. 0(500) GeV] weak gauginos, with ~ 1 — 2 TeV sfermion masses and weaker Higgs- 
stop couplings. In this case we find <5pmssm = 5.211 x 10~ 6 for the total payoff of the 
MSSM, whilst 5/9mssm = -1-338 x 10~ 5 for the specific SUSY Higgs part. As expected, 
the effects from the more massive SUSY degrees of freedom are relatively suppressed and 
the total contribution from sfermions and chargino-neutralinos is, in this particular case, 
comparable (albeit opposite in sign) to that of the Higgs bosons, so the overall MSSM yield 
is significantly smaller than in the previously discussed case. 

Let us now come back to the role of the leading higher-order corrections governed by 
enhanced Higgs boson self-couplings of the 2HDM. We spotlight significant effects which, in 
some particular corners of the parameter space, may well trigger relative deviations with 
respect to the one-loop prediction as large as A(5p eS )/5p^ ~ 0(30)% or above. Quite 
remarkably, these numerical results could be foreseen from the rough analytical estimate 

5o [1 - c{i] ~ 5o [1] ( 1 + 1 l Ahhh H 2 ( 4 7 ) 

r2HDM r2HDM \^ + 1 6 ^2 yp. J ' V*- ' > 



which follows from the approximate form factors quoted in Eq. (4.4). On the one hand, 
the triple Higgs boson self-couplings may reach maximum values of Ahhh ~ 0(1O 3 ) GeV, 
as we can read off Fig. [7] and Eq. fl2.18|) ; and, on the other, the typical Higgs boson masses 
considered in our analysis are of 0(100) GeV. Plugging these numbers into Eq. (4.7) we 
indeed retrieve A(5p eS )/5p^ ~ 50 - 100%, which agrees with the outcomes of our full 
numerical calculation to a remarkable extent. 

The parametric dependence on A5, as displayed in Figs. P-|l2|, exhibits the expected 
behavior bp ~ (-AA| + B A§ + C), which results from the structure of the (Born-improved) 
Higgs/gauge boson couplings gf^y ~ g (1 + 0(A§)) - entering of course square in the gauge 
boson self energies. It is for this reason that the leading correction to 5p is of 0(A|). 
Evidently this correction goes beyond the 2-loop level and therefore it tests in an effective 
way the largest possible effects that can be expected from the A5 parameter at all orders. 
This kind of effective approach is not new, let us note that it is fully within the line of 
reasoning presented in the work of Ref. [64], which was devoted to compute not only 
the full one-loop effects but also to test the largest possible quantum corrections that the 
trilinear Higgs boson self-couplings could produce at all orders in the pairwise production of 
neutral 2HDM Higgs bosons in a linear collider. As a matter of fact the one- loop diagrams 
used in that calculation are essentially the same ones that we are using to construct the 
effective Lagrangian interactions (|4.5| ) for the present higher order calculation, i.e. we are 
proceeding along the same philosophy as Ref. [64]. 

The leading large A5 corrections, if effectively realized in nature, would also manifest 
as a shift AMj eff] on the W-boson mass prediction. This shift we can roughly estimate 
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Figure 12: Analogous setup and results as in Fig. 0, for Sets V'-VIII' of Higgs boson masses with 
an inverted mass hierarchy [Af H ± — Mjsp < 0]. 
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Obviously this formula follows the same pattern as fl2.28|) up to the factor c^/s^ 



inherent to the definition of Ar^. The results for this estimate of AM, 



l-eff] 

w 



are displayed 



in the lower panels of Fig. ||. Notice that departures from the mere one-loop order by an 
amount of A{5p eS )/5p^ ~ O(30) % translate, for the mass sets under consideration, into 
a mass shift of up to A(5Mw) ~ MeV - as we can confirm from Table || and from 

the lower panels of Fig. |9[ The largest absolute higher-order corrections to 5p - attained of 
course for the maximum allowed (negative) values of A5 - are possible within non-SUSY- 
like Higgs boson mass spectra (Sets I and II). For the SUSY-like Sets III and IV, however, 
these effects are much more tamed, and amount to barely 0(1) MeV. Recall that despite 
the inconspicuous As-yield in this case, the full one-loop 2HDM correction on Ar is larger 
for those sets and therefore the overall impact is more sizeable (cf. third row of Table 2 
and Figs. 3-4). 

We also study the complementary Sets V-VIII and V'-VIII' (i.e. considering both the 



direct and the inverted mass hierarchies) in Figs. 11 and |12| respectively, and also in Ta- 
ble |||. For each of the mass sets we have different lower - and upper constraints on A5, that 
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-10.82 


0.93 


1.19 


20.22 
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-5.04 
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-6.41 


Set VII 


-8.04 
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0.72 


-3.37 


0.71 


-14.00 
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0.88 


1.27 


31.46 


12.07 


31.86 


12.38 


35.64 
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23.39 


5.82 


25.00 


2.57 


Set VIII' 


-5.40 


0.72 
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-34.04 


1.64 


-39.28 


1.68 


-40.62 


1.79 


-37.78 


1.77 



Table 4: Parameter space survey for Sets V-VIII and V'-VIIT in Table [jj The left-most columns 
quote the allowed A5 range, in agreement with the constraints from perturbative unitarity and 
vacuum stability. Alongside we indicate the corresponding value of tan /3 that maximizes the Higgs 
boson self- interaction enhancements for the largest allowed |Ag|. Within these optimal scenarios, 
and for different choices of a, in the right part of the table we quantify i) the relative size of the 
higher order effects on the <5p-parameter, A(Sp cS )/Sp^; and ii) the induced shift on the W-boson 
mass prediction [AM^ cff ], as defined by Eq. (p~8|). 

we explicitly account in the left-most columns of the table. The parameter tan f3 we fix 
accordingly, so that the Higgs self-coupling enhancements maximize for the largest allowed 
IA5I. The maximum attainable relative deviations on the <5p-parameter [A(5p eS )/5p^], 
alonside with the associated shift on the predicted W-boson mass [AM^ C ^], are docu- 
mented in the right columns. The results fall in the very same ballpark as for sets I- IV that 
we have examined previously. Again, the relative impact of the higher order effects can 
reach A(Sp eS )/5p^ ~ O(30%), but in this case dragging the W-boson mass predictions 
up to AM^ cff ' ~ 10 — 15 MeV in those regimes (viz. Sets V and V) which saturate 
the custodial symmetry bounds \5p\ < 10 -3 . For this reason, within these complementary 
mass sets we identify corrections to Mvv that can be up to three times bigger than in the 
original Sets I-IV. We should nonetheless point out that such remarkable effects beyond 
one-loop critically rely on a very large departure from the custodial symmetry limit, and 
hence are not favored by the EW precision data. Barring these borderline situations, a 
milder custodial symmetry breaking - e.g. at the level \Sp\ ~ 0(1O -4 ) - typically renders 
one loop corrections of up to AM-j^ eff ' ~ 40 MeV combined with higher order mass shifts 

[AM™- eff '] in the ballpark of 0(1 — 5) MeV. This result is nevertheless quite significant, as 
it implies that the characteristic higher order 2HDM corrections can be a few times larger 
than the pure bosonic 0(a1 w ) two- loop corrections (which perform at the < 1 MeV level), 
and even comparable to most higher order corrections within the SM. Finally, they can be 
as large as the estimated theoretical uncertainty (SM^ ~ 4 MeV) of the SM effects beyond 
the 2-loop order [38]. In all cases the results are seen to be very much responsive to A5, 
while they barely depend on the precise choice of the mixing angle a. 

The results illustrate once more the compromise between the two different trends that 
govern the overall impact of the higher order effects, and that somehow counterbalance 
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each other: i) heavier - versus lighter masses: the heavier the Higgs bosons, the weaker 
will be the loop-induced effects they generate - these shall be typically suppressed by inverse 
powers of such masses; ii) wider - versus narrower mass splittings; as we have discussed, 
these will determine the size of the pure one-loop piece [<5pW], and hence of the consequent 
mass shift [<5M^] and the corresponding higher-order contributions. This explains e.g. 
why Set V is the most responsive one, as it combines a sizable mass splitting between the 
CP-even, the CP-odd and the charged Higgs fields; and yet relatively light (viz. ~ 125 
GeV) neutral, CP-even states. 

5. Discussion and Conclusions 

In this work we have examined the traditional electroweak (EW) parameter Ar from the 
viewpoint of the general (unconstrained) Two-Higgs-Doublet Model (2HDM). That quan- 
tity, whose definition - confer Eq. ( |2.3| ) - was introduced more than thirty years ago 
[41,42]— spells out the quantum link between the EW gauge boson masses (M\y, Mz) and 
the Fermi constant (Gf), and thereby enables an accurate theoretical prediction for the 
mass of the W-boson [Mj^], based on the precise knowledge of the Z-boson mass [Mz], the 
fine structure constant [a] and Fermi's constant [Gf] as experimental inputs. Analyses of 
Ar have been instrumental in the past, both as precision tests of the SM and as a strategy 
to seek for hints of new physics, as well as to derive constraints on the parameter space of 
theories beyond the SM. In this paper, we have evaluated the novel 2HDM contributions 
(A^hdm) to the parameter Ar at one loop, namely the 0(a ew ) quantum effects that are 
determined by the non-standard Higgs boson corrections to the gauge boson self-energies. 
Let us note that to the best of our knowledge the full one-loop calculation of Ar within the 
2HDM was missing in the literature beyond the contribution to bp - the parameter which 
describes the breaking of the (approximate) SU(2) custodial symmetry. We have further 
folded the genuine full one-loop 2HDM effects on Ar with all the known SM corrections to 
that quantity. Of course the overlap of the SM Higgs boson contribution and the lightest 
neutral CP-even state of the 2HDM has been consistently removed. 

In the numerical analysis we have scrutinized a variety of regimes with phenomenolog- 
ical relevance and in full compliance with a brought-to-date set of restrictions that severely 
constrain the parameter space of the model - most significantly the unitarity and vacuum 
stability conditions, together with the low energy flavor physics observables B(b —> sj) 
and B° d - B d . Needless to say we have included as well the important phenomenological 
implications associated to the recent ~ 125 GeV Higgs- like state observed at the ATLAS 
and CMS detectors. We have examined various sets of Higgs boson masses, and used them 
to test the different response to considering: i) relatively light versus relatively heavy spec- 
tra; ii) relatively broad versus relatively narrow spectra; iii) a fully unconstrained choice 
versus a supersymmetric MSSM-like one. In all cases the genuine 2HDM effects on Ar 
exhibit a characteristic dependence on the Higgs boson masses, which manifests as a sort 
of antagonism between the following two tendencies: 1) the boosting of the 2HDM contri- 
butions to Ar from a significant mass splitting among the 2HDM Higgs states (caused by 
a stronger violation of the underlying SU(2) custodial symmetry); 2) the weakening of the 
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2HDM contributions to Ar when employing heavier sets of Higgs boson masses - obviously 
related to the suppression of the corresponding loop corrections. 

Let us also note that while our final results do of course depend on the details of 
the Higgs mass spectra, they are essentially unresponsive to the precise setup of Yukawa- 
coupling structures, i.e. they are insensitive to type-I or type-II 2HDM. This is not sur- 
prising at first sight, as no Higgs/fermion couplings enter the leading (one-loop) evaluation 
of the gauge boson self-energies. What is less obvious is that at higher orders, where the 
2HDM Yukawa couplings could, in principle, play a non-negligible role, they are finally im- 
potent, too, owing to the stringent theoretical and phenomenological constraints currently 
in force, in particular the condition that the regime of tan f3 = 0(1) is highly preferred. As 
a result the strength of the 2HDM Yukawa couplings is essentially the same as in the SM 
case within the physically allowed region. The different impact of type-I or type-II 2HDM 
on the evaluation of Ar enters our calculation in a completely indirect way, namely only 
through the specific (and significantly unequal) parameter space restrictions affecting each 
one of these models. 

Doubtlessly an expected but most salient feature of our analysis of Ar within the 
general 2HDM is the following: the sole addition of a second SUl(2) Higgs boson doublet 
significantly relaxes the persistent M$ — M^ p tension existing between theory and exper- 
iment. With only one Higgs doublet the global fits to EW precision data favor a Higgs 
boson lighter than expected (even below the LEP bounds) - precisely due to the depen- 
dence of Ar, and so of M$, on the Higgs boson mass(es). With an additional doublet 
the SM-like Higgs boson can stay comfortably at a mass value closer to the present exper- 
imental measurements since the rest of the (heavier) Higgs bosons provide the mechanism 
to compensate for the difference. All the explored scenarios feature very similar trends 
concerning both the obtained results and parameter space dependence for the quantities 
at stake. We identify a variety of regions within the 2HDM parameter space leading to 
typical variations of Ar in the ballpark of 5 (Ar) ~ O(10™ 3 ), which translate into W-mass 
shifts of 5M™ DM ~ 20 — 40 MeV. Being the current experimental error on the W-mass 
of 5M^ P = ±15 MeV (and the current discrepancy with the SM value of the same order 
of magnitude: M^ 1 - M^ p ~ -20 MeV), it is pretty obvious that the 2HDM effects are 
potentially very important to help the SM theoretical prediction (which is in deficit) to 
better match the experimental value. This is the more true if we recall that the future 
measurements may reduce the mass uncertainty further (hopefully at the ~ 5 MeV level), 
which means that we might be able to eventually confirm or exclude the SM value at a 
rather significant confidence level. 

While the above feature has been previously emphasized in the literature within the 
MSSM case, we confirm here for the first time that the general 2HDM shares this virtue. 
This may perhaps be viewed as not too surprising a posteriori since this fact basically 
depends (at one-loop) on the bulk contribution of the interactions of the gauge bosons 
with the new 2HDM Higgs states as compared to the single Higgs state in the SM. Let 
us note that in the context of the MSSM the Higgs sector is highly constrained and in 
general it provides a rather modest contribution to Ar. However, there we still have the 
additional one-loop contributions from the genuine super symmetric particles (e.g. from 
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squarks, sleptons, and chargino-neutralinos) . At one loop all these genuine SUSY effects 
are gauge-like interactions (since no Yukawa-like enhancement is possible at all at this 
order), and therefore they may amply compensate for the meager output from the MSSM 
Higgs sector, especially if the masses of some of the sfermions and/or chargino-neutralinos 
remain relatively light, i.e. at the few hundred GeV level. It is thanks to this feature 
that the overall theoretical predictions on Ar (and M\y) can be more in accordance with 
experiment within the MSSM than in the SM. 

What we have shown in this work is that this is also the case for the general non- 
super symmetric 2HDM, albeit for a completely different reason. Here the Higgs sector is 
less constrained and we have found that upon enforcing all the known restrictions from 
perturbative unitarity, vacuum stability, flavor and custodial symmetry, and direct search 
limits, the Higgs sector alone still possesses the capacity to effectively do the job, i.e. to 
render corresponding Ar and Myy predictions perfectly comparable to the total payoff of 
the MSSM. We believe it is useful and remarkable to have verified this fact explicitly, as 
it was not obvious a priori that the stringent current constraints on the 2HDM parameter 
space would still permit such possibility. 

No less remarkable in our study is the possible (and distinctive) quantum effects that 
the general 2HDM can provide beyond one-loop. This has been explored in the second part 
of our study, where we have dwelled upon the role of the higher order effects induced by the 
Higgs boson self-interactions on these EW precision quantities, specifically on the bearing 
they may have on the parameter Sp and hence on the characteristic Ar part carrying 
its influence (viz. —c^5p/s^). Owing to the various constraints, the Higgs boson self- 
interactions can only be enhanced in the limit of large values of the As-parameter in the 
general 2HDM Higgs potential. In turn the weak gauge boson self energies are responsive 
to them via Higgs boson- mediated corrections only at the two- loop level and beyond, 
which obviously makes the full task rather cumbersome. We have worked out a numerical 
estimate of such higher order effects by means of an effective Lagrangian approach which 
tracks the leading A5 self-coupling effects in the limit of large values of this parameter 
(within the allowed bounds). First, we have computed the dominant radiative corrections 
to all Higgs/gauge boson interactions in this limit. These corrections can be accounted 
for via a gauge invariant, UV-finite subset of (Higgs-mediated) one- loop diagrams from 
which we retain the 0(a ew A|) contributions. Second, we have rewritten these one-loop 
amplitudes as static form factors and used them to effectively dress the bare Higgs/gauge 
boson interactions. 

Quantitatively speaking, these higher order quantum effects on the 5p parameter may 
entail characteristic corrections to the 2HDM one-loop result as large as ~ 30 % - for 
Higgs boson self-interactions already bordering the unitarity limit. The corresponding 
impact on the theoretical prediction of the W-boson mass is typically of a few MeV. This 
is already a noteworthy result, since these typical effects are of the order of the present 
theoretical uncertainty within the SM prediction. Let us however mention that in limiting 
regions of the parameter space the higher-order 2HDM corrections could shift the W^-mass 
predictions up to ~ 10 — 15 MeV. Overall we find that the estimated higher-order effects 
within the 2HDM are small enough to remain beneath the limited precision of the current 
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EW data, but at the same time they could be at reach of the expected experimental 
precision attainable at the LHC and future TeV-range linear colliders. Finally, to better 
appreciate the potential significance of these effects, it is instructive to note that they can 
be not only substantially larger than the pure bosonic 0(a^ w ) two-loop contribution within 
the SM (which is long known to be of order of 1 MeV, at most) but they could even be of 
the same order than the remaining set of higher order pure QCD, EW and mixed QCD-EW 
effects. 

The upshot of our study is that the general 2HDM possesses, under the full set of 
currently known theoretical and phenomenological constraints, both the capability to im- 
prove the bulk theoretical prediction of the W-mass as compared to the SM, and also the 
possibility to provide distinctive quantum effects that could reveal its underlying dynamics 
through the role played by the 2HDM Higgs boson self-interactions. We have demonstrated 
that their characteristic impact on Ar might not be inconspicuous at all, and could in fact 
be quite relevant for an accurate theoretical determination of the EW precision quantities. 
Ultimately, they could emerge above the experimental precision and signal a smoking gun 
of physics beyond the SM and the MSSM 8 . Bearing in mind the remarkable improvement 
in the experimental accuracy which is expected for the forthcoming W-mass measurements 
at the LHC [viz. SM^ P ~ 10 MeV], this sort of analysis may become instrumental in the 
near future. Now that the LHC seems to be finally closing in on the Higgs issue, efforts 
on the theory side are of foremost importance to broaden the assortment of strategies that 
will be necessary to bring this research program to a most successful completion. 
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6. Appendix 

For the sake of completeness, we provide herewith a more detailed analytical account on 
selected aspects of our calculation. All UV divergences we handle by means of conven- 
tional dimensional regularization in the 't Hooft-Veltman scheme, setting the number of 
dimensions to d = 4 — e. As usual, we introduce an (arbitrary) mass scale \x in front of the 
loop integrals in order not to alter the dimension of the result in d dimensions with respect 
to d = 4. After renormalization (in the on-shell scheme, in our case) the results for the 
physical quantities are finite in the limit d — > 4. Furthermore, in the practical aspect of 
the calculation all one-loop structures are reduced in terms of standard Passarino-Veltman 
coefficients in the conventions of Ref . [94] . 

8 The potential importance of these distinctive quantum effects on Ar, as a trademark structure of 
(non-supersymmetric) extended Higgs sectors, was first suggested to the best of our knowledge in Ref. [64] . 
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One loop functions at zero momentum: the one-loop vacuum integrals that 
enter the evaluation of the parameter Sp, which is built upon the weak gauge boson 
self energies at vanishing momenta, cf. Eq, (|2.8|), read as follows: 
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where A e = 2/e + 1 — + log(47r ^ 2 ) and the function _F(x, y) is defined as follows: 



r£ ?-^Mf) x* v 



F(x,y) = F(y,x) = < 



x-y 



x = y 



(6.4) 



The tilde notation for the Passarino-Veltman functions, e.g. 



B(m 1 ,m 2 ) = B(0,m 1 ,m 2 ) , 



(6.5) 



indicates that these integrals are evaluated at zero momentum. The parameters A, B 
can be identified with the (squared of the) masses of the virtual particles propagating 
in the loop, A = m 2 , B = m|- 

With these expressions at hand, it is straightforward to write down a compact an- 
alytical form for 5p2HDM at one-loop in the 't Hooft-Feynman gauge, starting from 



the definition Eq. (gj): 
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—a 



cos 2 (P - a) F(M 2 , M 2 ± ) - F(M 2 , M 2 „ ) 



+ sin 2 (/?-a) F(M 2 o,M 2 ±)-F(M|o,M 2 ) + F(M 2 ,M 2 ± 



-3cos 2 (/3 - a) F(M 2 , + F{Ml„M 2 z ) - F{M^ ,M Z ) - F(M 2 , M&) 



(6.6) 



From the above equation we can explicitly read off how the size of bp depends on 
the mass splitting between the different Higgs bosons, as well as on the strength 
of the Higgs/gauge boson couplings - which is modulated by tan/3 and the mixing 
angle a. The first two lines of the full expression (|6.6j ) is the part that we have 
denoted <5p2HDM m ^ ec ' — > see Eq. ( 2. 15| ) . We remark that for M A o — > M H ± and 



|/3 — a\ — > 7r/2 (in which the h° field behaves SM-like) the full <5p2HDM — > 0. This is 
the precise formulation of the decoupling regime for the unconstrained 2HDM. 



In the case of the SM the Higgs contribution to the 5p-parameter fl2.8|) is not finite if 
taken in an isolated form. The complete bosonic contribution to Ar is of course finite 
and gauge invariant, and therefore unambiguous. To define a Higgs part of it is then 
a bit a matter of convention. What is important is that the complete M#-dependence 
is exhibited correctly and coincides in all conventions. After removing the UV-parts 
which cancel against other bosonic contributions one arrives at 



9 My , Mir 4 , , , , 

Mfln-f -M^ln-f -- (M 2 - M 2 ,) 



SM " 16 ^ 2 

+ F(M H ,M Z )-F(M H ,M W )} . (6.7) 

The explicit dependence on the scale p is unavoidable in quantities which are not UV- 
finite by themselves. It is however natural to set e.g. the EW scale choice p = Mw. 
In the limit > M^ we can see Veltman's screening theorem at work in the SM, 
as there remain no M 2 ter ms but a logarithmic Higgs mass dependence. Indeed, in 
that limit the expression ( |6.7D boils down to 



which coincides with the result quoted in Eq.( 2.12| ) of Sec 



The SM Higgs boson contribution to 5p can also be retrieved from the 2HDM result 
( |6.6D by selecting the h° parts of the contributions involving the h° and the gauge 
bosons, namely in the last line of that equation. By performing the identification 
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H = h° and removing the trigonometric factors we are led to 



dPsM ~ 16 ns 2 M 2 



F(Ml,M 2 z )-F(MlM, 



= n s- m — s- + ... (6.9 

16 vr 2 4 1 > 

We see that the last expression coincides with Eq. (|6.7| ) up to finite additive parts, 
which of course reflects the arbitrariness of the scale setting jj,. As we said, this is 
not important because the full bosonic contribution to Ar is finite and unambiguous. 
The fact that we can recover the SM result from Q6.6| ) in such a way suggests that the 
expression in the first line of ( p 7 ^ ) should be subtracted from ( |6.6| ) in order to compute 
the genuine 2HDM effects on 5p, i.e. the Higgs boson quantum effects beyond those 
associated to the Higgs sector of the SM. This is in fact the practical recipe that we 
follow in this paper. Finally, let us notice that the <5/C>2hdm P ar ^ °f (|6-6|) , i.e. the one 
which is completely unrelated to the SM Higgs contribution, is precisely the part of 
the full <5/?2HDM that violates the screening theorem in the 2HDM, as it is manifest 
from Eq. (fOp of Sec.[0. 



2HDM contributions to the gauge boson self-energies: We quote herewith 
their complete analytical form, in terms of the standard Passarino-Veltman coeffi- 
cients and following the conventions of Ref. [94]. The self-energies are evaluated 
for on-shell gauge bosons, e.g. p 2 = My [V = W , Z°], in the way they enter the 
calculation of Ar. 

— Two Higgs-boson contributions: 



| Higgs 
2HDM 



a 



-^ (M h 2 ) - A (Ml ) - A (Mlo) -2A (Ml, 



+4 cos 2 (/3 - a) Soo (M 2 ,, M 2 , M 2 ± ) + 4 sin 2 (/3 - a) B 00 (M 2 ,, M 2 , M 2 ±) 



+4 £oo (M^,M 2 ,M 2 ± ) 



(6.10) 



T 1 Higgs 
^ZJ 2HDM 



a 



16 7T S 2 ^ c\j 



-A (M? 



A (M? 



no) 



A (Ml ) - 2 (c^ - s^f A (M 2 ± ) 



+4 cos 2 (/3 - a) B 00 (Mf , M 2 , M 2 ) + 4 sin 2 (/3 - a) B 00 (Mf, M 2 , M 



A° ) 



+4 (c^ - s^Y Boo (M 2 Z , M 2 ± , M 2 d 



(6.11) 



Higgs/gauge boson and Higgs/Goldstone boson contributions 
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I Higgs/gauge 
2HDM 



|cos 2 (/3 — a) 
-cos 2 (/3 - a) Mw 



B 00 (M& , M 2 o , ) - B 00 (Afw » , ) 



S (M^ , M| , ) - Bo (Mw , M 2 , ) 



}• 

(6.12) 



I 1 Higgs/gauge 
jZ J 2HDM 



{cos 2 (/3 - a) \Boo [Ml, M 2 „, M|) - £ o (M|, M 2 , M 2 ) 
cos 2 (/3 - a) Ml [b (M|, M 2 , M|) - B (Mf, M h 2 D , Mf )] } . (6.13) 



Let us notice that, in the last two expressions, we have explicitly removed the overlap 
with the SM Higgs boson contribution, to wit: 



2HDM 



= E V - Xy] SM « 

cos\f3-a)Boo (M V ,M { 
-M^ cos 2 (fi - a)B( 
- B { " • • 



Boo (M 2 ,M 2 ,M 2 ; 



r 2 o,M 2 ) + sin 2 (/3- 

? (My, Mlo,My) - My S hl 2 (/? - a) B ( 

| Hsh0 + My Bo (My, M 2 , M 2 )| H=h() 



a) 5, 

2 



? 00 (My, Ml,, My 



?o (M 2 ,M 2 ,M 2 ) 



(6.14) 



Effective Higgs/gauge boson interactions To better illustrate how we build up 
in practice the effective Higgs/gauge boson couplings employed in this study, here- 
with we provide explicit analytical details for the construction of one of such Born- 
improved interactions. The calculation we carry out with the help of the standard 
algebraic packages Feyn Arts and FormCalc [89, 94] . Without loss of generality, 
let us take the concrete case of the Z boson coupling to the CP-odd and the light CP- 
even neutral Higgs bosons [ShOA°z°]- A sample of the Feynman diagrams describing 
the 0(A 2 ) corrections to this coupling is displayed in the upper row of Fig. [j| The 
general structure of the associated form factor a h o A o z o may be casted as: 



a h o A o z o 



ReVL. 



1 



1. 



tan(/3 



p 2 =0 



A° 



wave— function 



M 2 



h°H0 



mixing 



(6.15) 



Notice that we define our form factors to be real, in order to preserve the hermiticity 
of the effective Lagrangians derived from them. The different building blocks of 
Eq. ( |6 .15 ) correspond to: 

a) V^oaoz ) the genuine vertex corrections (cf. e.g. the first two diagrams in the upper 
row of Fig. |8|). For illustration purposes, we provide its complete analytical form: 
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V^o A o z o (0) = |A H o A OAoA h o h o H oC'2(M^ ,M^o,M^o) + 

A h o A o A oA h o A o A o Ci(M^p, Mjp, M%p) - A h o A o A oA h o h o h o C 2 (M A o, M 2 , M 2 ) 

+ tan (/J - a) A h o A o A oA h o h o H o C\ (M A0 , M 2 , M^ ) 

— Ah°a°a° AhOA«A° Ci(Mj i p, M^p, M^ ) 

+ Ah°A°A AhOH°H° C2 (M A o , M^o , M^o ) I • 

(6.16) 

b) the wave-function corrections associated to each of the external Higgs boson legs 
(including, as we single out in the last term of Eq. ( |6.15| ) , the h° — H° mixing one-loop 
diagrams): 

Sho(o) = {awhob;(4.4) + a2 0h±h± ^(m2 ± ,m^ ± )+ 

AhO H o H o B' (M^ , M^o) + Aho A o A o £?o(M|p, M A o) + A^o h o h o B' {Mlo, M 2 ) j . 

(6.17) 



1 

+ 2L 



S Ao(°) = _ Y6^2 ^hOAOAO^o( M AO' M ho) + ^HOAOAO-^o( M HO) M Ao) • ( 6 - 18 ) 



S h o H o(0) = ^ { A h"ao A o A h o A o A o [5 (M|„,Ml„) - KeB (g 2 , M 2 , M 2 )] 
+A h o h o h o A h o h o H o \B {M 2 ,M 2 ) - 5fte B (q 2 , M 2 , M 2 )] 
+A H o H o H o A h o H o H o [S o (0, M 2 , M 2 ) - 5 (g 2 , M 2 , M 2 )] 
+2 (A h o H o H o A h o h o H o [S (M 2 o, M 2 ) - 3Je5 (g 2 , M 2 „ M 2 )] 
+A h o H+H - A H o H+H - [So(M 2 ± ,M| ± )-5R eJ Bo( 9 2 ,M 2 ± ,M 2 ± )])}. (6.19) 

In the last equation, the h°— H° mixing self-energy S h OHo((? 2 ) = S h ono((? 2 )+(5 Z h ono(q 2 — 
M? )/2 + <5 Z h OHo((? 2 — M^ )/2 — 5M^ 0h0 involves the renormalization of the mixing 

angle a, that we an chore via the relation 3f?e£h H (<? 2 ) = according to [64], with 
the renormalization scale chosen at the average mass q 2 = (M^ + M^ )/2. As 
mentioned above, the tilded Passarino-Veltman functions are evaluated at vanishing 
external momentum. 

Let us note in passing that, for the case of the ffhw-type couplings, and due to he 
fact that just one single scalar leg is present there, only pieces of type b) shall give 
rise to 0(X 2 ) contributions. The same holds as well for the Higgs/gauge/Goldstone 
boson couplings [ghvc]- 
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